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Abstract. 

We review recent results concerning the renormalization group (RG) transformation 
of Dyson's hierarchical model (HM). This model can be seen as an approximation of a 
scalar field theory on a lattice. We introduce the HM and show that its large group of 
symmetry simplifies drastically the blockspinning procedure. Several equivalent forms 
of the recursion formula are presented with unified notations. Rigorous and numerical 
results concerning the recursion formula are summarized. It is pointed out that the 
recursion formula of the HM is inequivalent to both Wilson's approximate recursion 
formula and Polchinski's equation in the local potential approximation (despite the 
very small difference with the exponents of the latter) . We draw a comparison between 
the RG of the HM and other RG equations in the local potential approximation. The 
construction of the linear and nonlinear scaling variables is discussed in an operational 
way. We describe the calculation of non-universal critical amplitudes in terms of 
the scaling variables of two fixed points. This question appears as a problem of 
interpolation between these fixed points. Universal amplitude ratios are calculated. 
We discuss the large- limit and the complex singularities of the critical potential 
calculable in this limit. The interpolation between the HM and more conventional 
lattice models is presented as a symmetry breaking problem. We briefly introduce 
models with an approximate supersymmetry. One important goal of this review article 
is to present a configuration space counterpart, suitable for lattice formulations, of 
other RG equations formulated in momentum space (sometimes abbreviated as ERGE). 
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RG: Renormalization Group 

ERGE: Exact Renormalization Group Equations 

LPA: Local Potential Approximation 

HM: (Dyson's) Hierarchical Model 

HT: High- Temperature 

UV: Ultra- Violet 

IR: Infra-Red 

£: the change in linear scale after one RG transformation 

the sum of all the fields 
Q: the total number of sites 
(pc- the classical field or magnetization density 
0„: the sum of the fields in a block of size 2" 
Wn the "raw" local measure after blockspinning 2'^ sites 
Wn'- the rcscalcd local measure after n RG transformations 
Tin'- Baker's form of Wn 
Qn- the Gallavotti's form of >V„ 
c = 21-2/^ 
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1. Introduction 

Quantum field theorists face tlie arduous task of figuring out tlie large scale implications 
of models defined by interactions at a small scale. In general, the large distance 
behavior of the theory can be encoded in an effective action Sgff which describes the 
interactions of the zero-momentum modes of the fields. Knowing Seff, we can answer 
important questions regarding, for instance, the stability and triviality of the theory. 
The renormalization group (RG) method [1, 2, 3, 4] was designed to calculate Seff by 
a sequence of small steps where the high energy modes are integrated progressively. 
This procedure generates a sequence of Seff, a where A is the scale above which we have 
integrated all the modes. In the process, many new interactions are introduced and with 
them functions also evolving with A. This can be seen as a flow in the space of theories. 
We call this type of flows the RG flows. Important simpliflcations may occur near flxed 
points that have only a few unstable directions and universal properties may emerge. 
Despite the conceptual beauty of the construction, the practical calculation of the RG 
flows remains very difficult. The behavior near flxed points can usually be handled by 
some expansion (weak coupling, strong coupling, etc.), but interpolating between two 
flxed points is in general a major difficulty. Unfortunately, this is essential to calculate 

^eff- 

In the following, we review our present understanding of the interpolation between 
flxed points for Dyson's hierarchical model [5, 6], a model for which the calculation 
of Seff, A can be reduced to the calculation of the effective potential Ve//,A and can 
performed numerically with great accuracy. Other type of hierarchical models have been 
discussed in the literature (see, for instance, reference [7]), however, in the following, 
we only consider Dyson's model. The Hierarchical Model (HM) will only refer to 
Dyson's model in this review. The RG transformation for the HM can be expressed 
as a simple one-variable integral equation very similar to the so-called approximate 
recursion formula proposed by Wilson [3]. In many respects, the techniques involved 
in the solution of this model can be compared to those used in elementary quantum 
mechanics. The reason for this remarkable simplicity is that the kinetic term of the 
model is not renormalized. In other words, it is a model for which the Local Potential 
Approximation (LPA) is exact. 

In recent years, the LPA has been widely used in the context of Exact 
Renormalization Group Equations (ERGE) and has generated a lot of interest. We 
recommend references [8, 9, 10] for reviews of the recent progress. ERGE allow in 
principle the study of global and nonlinear aspects of the RG flows of fleld theoretical 
or statistical models. However, truncation methods are necessary in order to make 
practical calculations. One particularly popular choice is to combine Polchinski's ERGE 
[11] with the LPA [12]. This results in a simple partial differential equation (called 
"Polchinski's equation" below) for the effective potential. Polchinski's equation can also 
be obtained [13] as an inflnitesimal version of the RG equation for the HM. This suggests 
that the linearized theories should be close to each other. Accurate calculations of the 
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critical exponents [14, 15, 16, 17, 18, 19] show that the exponents differ only in the 
fifth significant digit. In reference [18], it was believed that the two exponents should 
coincide, however this is not the case (this will be explained in section 9). 

Originally, the HM was introduced by Dyson [5] as a long-range ferromagnetic 
Ising model with couplings weaker than the one- dimensional Ising model with long- 
range couplings falling off like a power. The existence of a phase transition and of 
the infinite volume limit for some range of the parameter controlling the decay of the 
interactions, can be proved rigorously. Historical motivations and rigorous results are 
discussed in more detail in section 6. The model was rediscovered a few years later by 
Baker [6] in an attempt to construct models for which Wilson's approximate recursion 
formula or an integral equation of the same form become the exact RG transformation. 
In this context, the HM appears rather like an approximation of a scalar field theory on 
a D-dimensional lattice. Later, we call the process of interpolation between the HM and 
lattice models with nearest neighbor interaction the "improvement of the hierarchical 
approximation" . 

One important goal of the review article is to stress the similarities between 
Polchinski equation (a generic computational tool) and the HM (often perceived as 
a toy model). The motivations for using the HM model are: 

• It is a lattice model right from the beginning. 

• High-accuracy methods exist to calculate numerically the critical exponents and 
the RG fiows for arbitrary potentials. 

• Conventional expansions (weak and strong coupling, e-expansion) can be 
implemented easily and large order series can be obtained. 

• The hierarchical approximation can in principle be improved [20, 21]. 

The study of critical phenomena has reached a stage where many methods have been 
refined to a point where they provide numbers very close to each other [22, 23]. In 
the case of the HM, all the approximations can be compared to very accurate numerical 
answers. The ability to construct the RG fiows very accurately means that we can study 
general features of these flows far away from fixed points. This type of study is also 
possible using the LPA [24, 25, 8, 26]. 

An interesting feature of the HM is its discrete scale invariancc which, depending 
on the context, can be seen either as an annoyance [27, 28] or an interesting intrinsic 
property [29]. The review is focused on doing calculations directly in 3 dimensions. The 
e-expansion near 4 dimensions, which is reviewed in reference [30], is not discussed here. 

The paper is organized as follows. In section 2, we briefly review the basic 
concepts of the RG method, the scahng hypothesis and the practical difficulties of 
blockspinning. This section motivates approximations that simplify the blockspinning 
procedure. Dyson's HM is introduced in section 3, with the notation used later in 
the review. Several equivalent forms of the recursion formula are presented in section 4. 
More general recursion formulas which coincide with the HM's one for a particular choice 
of what we call the scale parameter, are introduced in section 5. This includes Wilson's 
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approximate recursion formula [3] which is inequivalent to the recursion formula of the 
HM. In the following, we denote the scale parameter £. In the literature, the same 
parameter is often denoted L, but we preferred to keep that symbol for the linear size 
of the whole system. 

Having introduced the basic concepts, motivations for the HM, numerical and 
rigorous results are reviewed in section 6. The numerical treatment of the recursion 
formula is discussed in section 7 and perturbation theory with a large field cutoff in 
section 8. The relation to Polchinski's ERGE in the LPA is discussed in section 9. 

We introduce the non-linear scaling variables associated with a fixed point in section 
10. These quantities, originally introduced by Wegner [31], transform multiplicatively 
under a RG transformation. They have features similar to the action-angle variables 
used in classical mechanics. Small-denominator problems are in general present and 
need to be discussed in each case. In section 11, we describe the non-universal critical 
amplitudes as RG invariants made out of two non-linear scaling fields. A practical 
calculation based on this method requires the ability to use both sets of scaling fields in 
a common intermediate region, in other words, to interpolate between the fixed points. 

The notion of a nontrivial continuum limit, originally introduced by Wilson [32] 
is reviewed and applied to the HM in section 12. Calculations of critical amplitudes 
and their universal ratios are then discussed. In section 13, we introduce the extension 
of the HM for N components and discuss the large- limit. We compare the results 
with those obtained with Polchinski equation emphasizing the difference in the 5th digit 
already mentioned above, which reflects the non-equivalence of the models. 

The possibility of improving the hierarchical approximation by breaking its 
symmetries in a systematic way is discussed in section 14. Compared to the 
improvement of the LPA by the derivative expansion, the improvement of the 
hierarchical approximation is an underdeveloped subject. On the other hand, it is clear 
that much progress remains to be done in the ERGE approach in order to match the 
accuracy of other methods [33, 34, 23] for the calculation of the critical exponents [35]. 
We hope that this review will facilitate the communication between the two approaches. 
Reference [36] is a very recent example of progress made in this direction. Finally, 
the possibility of including fermions in approximately supersymmetric models is briefly 
discussed in section 15. 

One motivation to study global and nonlinear aspects of RG flows not covered in 
this review is to improve our understanding of quantum chromodynamics. In this theory 
of strongly interacting quarks and gluons, weakly interacting particles are seen at short 
distance (asymptotic freedom), while nonperturbative effects cause the conflnement 
of quarks and gluons at large distance. Understanding how these two behaviors can 
be smoothly connected is a single theory amounts to interpolate between two flxed 
points of the Renormalization Group (RG) transformation. Despite recent progress 
[37, 38, 39, 40, 10, 41, 42] it remains a challenge to understand conflnement in terms of 
weak coupling variables. This is work for the future. 
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2. General and practical aspects of the RG method 

In this section, we present the basic ideas behind the RG method and point out the 
practical difficulties associated with the so-called blockspin method. 

2.1. Statement of the problem 

We are interested in the large distance (low momentum, long wavelength) behavior of 
scalar models. We consider a generic scalar model with a lattice regularization and a 
action S. The scalar field (f)x is coupled linearly to a constant source J. We call the 
total field 

* = E'^-' (2.1) 

X 

and the total number of sites Jl. With these notations, the partition function reads 

/+00 r+oo 
... J](i0,exp(-5 + $J) , (2.2) 
oo •/— oo J. 

We define the average at J = of an arbitrary function A of the fields as 

/+0O r+oo 
... lldcl>,AeM-S)/Z[0] (2.3) 
■oo J —oo ^ 

It is clear that 

oo 

Z[J]/Z[0] = 1 + 5^ - < > (2.4) 

g=l ^' 

The connected parts can be obtained by taking the logarithm of this expression. At the 
lowest orders, 

< $2 >c = < $2 > -(< $ >)2 

< $3 = < $3 > -3 < $2 >< $ > +2(< $ >)^ 

< $4 ->c ^ < $4 > -4 < $3 >< $ > -3(< $2 >f 

+ 12 < $2 > (< $ >)2_g(^< $ >)4 (^2.5) 

In general, we expect that 

< >"(x Q , (2.6) 

and we define 

^(g) >c (2.7) 

even though the individual terms scale faster than Q. Unless wc take the infinite volume 
limit, we should in principle write x'n order to remind the dependence on the volume. 
In the rest of this subsection, this dependence will be kept implicit. We then write the 
generating function of the connected densities 

oo ^ 

{lM\n{Z[J]/Z[0]) = J2-^J'x^'^ (2-8) 
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It is common to call 

G[J] = -(l/Q)ln(Z[J]/Z[0]) , (2.9) 
the Gibbs potential and it is clear that 

X^i^ = -d'iG[J]/dJ'i. (2.10) 
We now introduce 1 in the functional integral in the following way 

/ + 00 

d(t>,8{(t>c-^l9) . (2.11) 
-oo 

The partition function becomes [43] 

Z[J]^ p d0eexp(^-O(K//(0c)-0cj)^ , (2.12) 



with 

f+OO 



/+00 r+oo 
... J]d0,exp(-5)5(0e-$/O) . (2.13) 
-oo J —oo J. 

If for any J, T4y^j(0c) — <l>cJ has a unique minimum at 0c — <Pc, we have in the limit of 
arbitrarily large that 

Z[J] oc exp - n(Veffi^c) - ^cj)^ . (2.14) 

0c is a function of J defined implicitly by 

9K//(0c)/90cU.=<^-. = J . (2.15) 
G[J] and Veff{(l)c) are related by a Legendre transform 

G[J]=Veff{^c)-^cJ (2.16) 

In the ferromagnetic language, (pc is called the magnetization and J the magnetic field. 
In analogy with the gas-liquid transition, V^ff plays the role of the Helmholtz potential, 
0c the role of minus the volume and J the role of the pressure. When we are in the 
gas phase below the critical temperature (but above the triple point temperature) , if we 
increase the pressure keeping the temperature fixed, the volume decreases until a critical 
pressure is reached where the gas and the liquid can coexist at an equilibrium pressure 
but with different specific volumes. The discontinuity in the magnetization is analog to 
the specific volume discontinuity. The sign is justified by the fact that as we increase 
J, the magnetization increases, but if we increase the pressure, the volume decreases. 

In the simple case where Vg/f is an even function with a unique minimum at zero, 
we can expand 

oo ^ 
9=1 

Using the derivative of this expansion to express J in terms of 0c, plugging into equation 
(2.16), using equation (2.8) and solving order by order in 0c, we obtain the well-known 
relations 

r(2) = 1/^(2) , 
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r(4) ^ _ ^(4)^(^^(2)^)4 ^ 

r(6)^10(;^(^))V(;^(2))V_^(6)/(^(2))6^ (2.18) 

and so on. The calculation of the 1^// is an important objective the RG method. 
2.2. Basic aspects of the RG method 

We now consider a scalar model on a D-dimensional cubic lattice with a lattice spacing 
a. The RG transformation proceeds in two steps. First, we integrate the £^ fields in 
blocks of linear size £a while keeping the sum of the fields in the block constant. We 
then divide the sums of the fields by a factor £(2+^-»^/)/2 Q^^id treat them as our new 
field variables. The exponent rj is introduced in order to keep the canonical form of the 
kinetic term and its calculation is non-trivial. We then obtain a new theory in terms of 
a new field variable which is equivalent to the previous one as long as we only consider 
processes involving energies smaller than the new ultraviolet cutoff ~ 1/ {£a) . Given an 
original action, we assume that the RG transformation provides a new effective action 
expressed in terms of the new field variable. We postpone the discussion of the practical 
aspects of the partial integration to subsection 2.3. 

The information that is kept during the RG transformation is encoded in the average 
values of all the integer powers of the sum of the fields in the blocks. We call these 
average values the "zero momentum Green's functions at finite volume". This set of 
values can be thought of as an element of an infinite vector space. In the following 
we call the trajectories in this space the RG fiows. If we start with an infinite volume 
lattice, the RG transformation can in principle be repeated an infinite number of times. 
If this can be done, the resulting effective theory is limited to the zero-momentum 
Green's functions. We are particularly interested in finding the fixed points of the RG 
transformation, since the above limit can be simplified in the neighborhood of a fixed 
point. 

The RG transformation has several obvious fixed points. If the interactions 
arc limited to quadratic ones and there is no restriction on the range of the fields, 
the partition function can be calculated exactly using diagonalization and Gaussian 
integration. The model describes non-interacting particles of a given mass. We call this 
mass in cutoff units m. It can be determined in cutoff units using the zero-momentum 
two point function. 

After one RG transformation, the mass for the new effective theory is mi. The 
Gaussian fixed point (also called the trivial fixed point) corresponds to the theory with 
m — 0. There are usually other fixed points which corresponds to stable phases. For 
instance, the high-temperature fixed point which can be thought as an infinitely massive 
theory where the fiuctuations about the zero field value are entirely suppressed. 

In addition to these obvious fixed points, we expect non-trivial fixed points when 
D < 4. These are characterized by one or more unstable direction. The critical 
hypersurface is given as the stable manifold (e.g. the basin of attraction) of this non- 
trivial fixed point. Its codimension is the number of unstable directions. In the rest of 
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this section, wc specialize the discussion to the case where there is only one unstable 
direction. The stable manifold can then be reached by considering a family of models 
indexed by a parameter which can be tuned in order to cross the stable manifold. In 
field theory context, we usually pick the bare mass to accomplish this purpose. In the 
statistical mechanics formulation, the inverse temperature /? can be tuned to its critical 
value which is a function of the other interactions. 

Near the non-trivial fixed point, we can use the eigenvectors of the linearized 
RG transformation as a basis. As far as we are close to the fixed point, the average 
values of the powers of the rescaled total field stay approximately unchanged after one 
transformation. However at each iteration, the components in the eigendirections are 
multiplied by the corresponding eigenvalue. In the following, we denote the eigenvalue 
larger than 1 as Ai. As we assumed that there is only one unstable direction, there is 
only one such eigenvalue and it is the largest one. 

After repeating the renormalization group transformation n times, we have replaced 
sites by one site and associated a block variable with it. We define the finite 
volume susceptibility x^^^ as the average value of the square of the sum of all the 
(unrescaled) fields inside the block divided by the number of sites We can estimate 
X^^^ near the non-trivial fixed point. The average of the square of the rescaled variables 
is approximately a constant that we call Ki. To get the susceptibility, we need to go 
back to the original variables (so we multiply Ki by £(2+£>-r))n^ g^^j divide by the volume 
£Dn^ We can also take into account the motion along the unstable direction in the linear 
approximation. In summary, near the fixed point, 

ri2-v)^K^ + - . (2.19) 



(2) 
n 



The constant K2 depends on the way the critical hypersurface is approached when 
f3 is varied close to f3c- equation (2.19) is valid only if the linearization procedure is 
applicable, in other words if A"(/3c — /?) << 1- On the other hand, when n reaches some 
critical value n* such that 

Xf {(3,-/3) c^l, (2.20) 

non-linear effects become important and the sign of {(3, — (3) becomes important. In the 
following we consider the case of the symmetric phase {j3 < j3c) which is simpler. The 
order of magnitude of x starts stabilizing when n gets of the order of n* and 

X^^«r*M«(/?c-/?)-\ (2.21) 



with 



and 



7 =(2 -77)^, (2.22) 



u = ln£/lnAi . (2.23) 

When n >> n*, the trajectories fall into the completely attractive HT fixed point. 
The approach of the fixed point is characterized by corrections proportional to negative 
powers of the volume. 
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Wc can also take into account the corrections due to the so-called irrelevant 
directions which correspond to eigenvalues less than 1 and eigenvectors along the stable 
manifold. We call the largest of the irrelevant eigenvalues A2 < 1. As n increases, the 
distance along the corresponding eigenvector shrinks as Ag. The relative size of the 
correction should be proportional to 

Af ^ - (3)^^ , (2.24) 

with 

A, = -InAa/lnAi . (2.25) 

As has a lower index s which is short for subleading and should not be confused with 
the gap exponent that will be denoted A^. In summary, in the large volume limit and 

for /3 ^ /3~ , we have the parametrization 

X^^) ~ (/3c - Py^Ao + Ai(/?c - + ...]. (2.26) 

Linearized calculations near the non-trivial fixed point allow us to calculate the 
exponents 7 and A,^, but not the amplitudes Aq and Ai. 

It should be emphasized that often, the lattice spacing docs not appear explicitly 
but is considered as a function of the bare parameters. The relation between the two is 
usually established by relating lengths calculated in lattice spacing to physical lengths. 
The continuum limit consist in taking a trajectory in the space of bare parameters which 
corresponds to the limit of zero lattice spacing. 

When P > Pc, we are in the low temperature phase and x^^^ 7^ 0. The sign of 
X^^^ 7^ is positive (negative) if we take the hmit J — > by positive (negative) values. 
Additional subtractions are then required as shown in equation (2.5). The case of the 
low temperature will be discussed in subsection 7.5 for the HM. 

Relations among exponents can be obtained from the scaling hypothesis [44, 45, 46]. 
Making explicit the dependence of the non-analytical part of the Gibbs potential, 
denoted Gg, on the reduced temperature t = {Pc/ P) ~ 1; the scaling hypothesis amounts 
to have 

Gs[e^'t, i^'^J] = i^Gs[t, J] . (2.27) 
The exponents y are connected to the critical exponents discussed above by the relations 



V 



(2.28) 

£' + 2-77 



2 

It is also common to introduce exponents which describe the correlations as a function 
of the inverse distance for the conjugated variables at criticality, namely 

x^, = D-y,, (2.29) 

Xt=D-yt, (2.30) 

for the field and for the energy respectively. The scaling hypothesis will be further 
discussed in subsections 3.3, 7.4 and 7.5. 
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2.3. Practical aspects of blockspinning 

In the previous subsection, we have assumed that it was possible to integrate over the 
variables in a block while keeping the sum of the field constant. This procedure is 
usually called "blockspinning" , an idea which can be traced back to Kadanoff [47] . In 
practice, blockspinning is usually quite complicated. However, for models with actions 
quadratic in the fields, it is possible to do it analytically. Of course, for such models, all 
the correlation functions can be calculated exactly, but the procedure can be considered 
as the first step in a perturbative expansion. To fix the ideas, we can consider a one- 
dimensional lattice model where the Fourier transform of the two point function (the 
propagator) is G{k). Assuming that the lattice sites are labeled by integers, we have the 
periodicity G{k + 27r) = G{k). If we now partition the lattice into blocks of even-odd 
pairs of neighbor sites and blockspin within these blocks, we obtain a new two-point 
function on a new lattice with a lattice spacing twice larger. If we denote the Fourier 
transform of the two-point function after n steps Gn{k), the iteration formula is: 

Gn+iik) = (1 + cos{k/2))Gn{k/2) + (1 - cos(A;/2))G„(A;/2 + tt) . (2.31) 

It is clear that the 27r periodicity is preserved and no sharp edges are introduced at 
least in a finite number of iterations. The construction can be extended in arbitrary 
dimensions and was used as the starting point for the finite- lattice approximation [48] . 

In order to get an idea of the difficulty to extend the procedure when higher order 
interactions arc introduced, it is instructive to consider the simple case of a 4-sites ring 
with nearest-neighbor quadratic interactions and local quartic interactions. 

/ + 0O /■ + 0O 

. . . / d(f)id(l)2d(j)3d(f)^e-^^-^^-^' , (2.32) 
-OO J —OO 

with nearest neighbor quadratic terms 

= (01 - 02)' + (03 - 04)' , (2.33) 

and 

= (02 - 03)' + (04 - 0l)' , (2.34) 
and local quartic interactions 

^4 = 0t + 0^ + 0^ + 0l . (2.35) 

For further convenience, we also define the next to nearest neighbor quadratic 
interactions 

S^^^ = (01 - 03)' + (02 - 04)' (2.36) 

We now try to blockspin. We pick (1,2) and (3,4) as our basic blocks and then 
combine them into the block containing all the sites. More specifically, it consists in 
introducing 1 in the integral in the following form: 

/+00 r+oo n+oo 

I I d$d$(i,2)(i$(3,4)(^(^ - ^(1,2) - ^(3,4)) 
■OO J —OO J —OO 

X ^($(1,2) - 01 - 02)5(*(3,4) - 03 - 04) (2.37) 
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If we can perform the integration over 0i — 02 and 03 — 04, and then over $i — $2, we 
will be able to write 

/+00 
d^e-^^ff-^^^ . (2.38) 
■00 

In this illustrative model, SeffX^) is the main quantity of interest. Because of the 
quadratic terms S*^, we cannot perform the integral over 0i — 02 independently of the 
integral over 03 — 04 and blockspinning is potentially more difficult than evaluating 
the integral without intermediate steps. It is interesting to notice that the terms of S2 
connect the fields across the blocks and are not invariant under independent interchanges 
(1 <-> 2) and (3 <-> 4). These two transformations are of order 2 and commute, they 
generate a group of order 4. If we average S2 over the 4 elements g of this group, we 
obtain 

(1/4) (y) = (l/2)(5f + ^2™) (2.39) 

9 

If we replace S2 by this average, we have twice more terms but with half of the strength. 
Also, the new terms have a longer range. At first sight, it is not clear that the situation 
is better than before. However, if we combine this average with one half of the harmless 
S2, which is invariant under the above mentioned symmetry group, we obtain 

(l/2)(52^ + 5f + ^2''^^) =4(0? + 0^ + 0^ + 02)-(0, + 02 + 03 + 04)2 .(2.40) 

The first term affects only the local measure and the second can be incorporated directly 
into SeffX^)- This simple example illustrates how a symmetry can be used as a guide 
to build a (modified) model where blockspinning is feasible. 

Dyson's hierarchical model is a model where a hierarchical exchange symmetry 
among the sites is built-in and allows to blockspin the partition function by performing 
a sequence of one-dimensional integrals. 

3. Dyson's hierarchical model 
3.1. The Model 

In this subsection, we describe Dyson's hierarchical model with the notations used in 
most of the rest of this review. The relationship between this formulation and other ones 
found in the literature is discussed in sections 4 and 5. The model requires 2""""^ sites. 
We label the sites with rimax indices x„„^^^, xi, each index being or 1. Wc divide the 
2"'"""= sites into two blocks, each containing 2"'™-''="~^ sites. If a:„^^^ = 0, the site is in the 
first block, if Xn^ax = the site is in the second block. Repeating this procedure n times 
(for the two blocks, their respective two sub-blocks , etc.), we obtain an unambiguous 
labeling for each of the sites. The indices on the left provide the coarser division while 
the indices on the right provide a finer division. With an appropriate choice of origin, 
the indices can be interpreted as the binary representation of the site numbers. This is 
represented graphically in figure 1: 
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The nonlocal part of the total energy reads 

^ = -2E(|r E ( E '^(-.a..--))' • (3-1) 

The index n, referred to as the "level of interaction" hereafter, corresponds to the 
interaction of the total field in blocks of size 2". The constant c is a free parameter 
assumed positive and which controls the decay of the iterations with the size of the 
blocks and can be adjusted in order to mimic a D-dimensional model. This point is 
discussed in more detail below (sec cquation(3.23)). 

The field is integrated over a local measure which we need to specify. 

In the following, we will often work with the Ising measure, VFo(0) = ^(0^ — 1) or a 
Landau- Ginsburg measure of the form VFo(0) = exp(— im^0^ — A0"^). 

In the case of the Ising measure, the only free parameter is c. If all the (p = +1, 
then the cost in energy for flipping one spin is 2 ^;^^(c/4)'(2' — 1) and is finite in the 
infinite volume limit only if c < 2. In the following it will be assumed that < c < 2. 

The hierarchical structure of equation (3.1), allows us to integrate iteratively the 
fields while keeping their sums in blocks with 2 sites constant. 

W^n+i(0(i,2)) = exp ((/5/2)(c/4)"+V?i,2)) X J^d^Wni^+OWni^-O ,(3-2) 

where 0(i,2) is understood as the sum of the two fields in the block. After rimax 
integrations, we obtain 

^^W^n_(^)=e-^^«//(*/"), (3.3) 

with 

Q = 2'*'""- (3.4) 
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Figure 1. Divsion in bloclcs for 8 sites. Tlie lower row is X3X2X1. 
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Remarkably, the symmetries of the model have allowed us to calculate the partition 
function and the effective potential by calculating only rij^ax independent integrals 
instead of the 2"™""= coupled integrals that one would naively expect. 

The symmetry group of the HM is of order 2^^^, which is one half of the number 
of configurations for a Ising measure. This can be seen as follows. In any of the blocks 
of size 2 we can interchange two sites. The exchange can be done independently in each 
of the blocks and this symmetry can be seen as local. There are Q/2 blocks of size 2 
and so 2^/^ distinct symmetry transformations. Similarly, we can interchange the blocks 
of size 2 inside any block of size 4 which generates 2^/^ independent transformations. 
We can similarly generate independent symmetries until we reach the group of order 2 
exchanging two blocks of size fl/2. Using 

1/2 + 1/4 + ... + (1/2)"-"^ = 1 - 1/Q , 

we obtain the announced result. 



3.2. The RG transformation 

As explained in subsection 2.2, the RG transformation consists in a blockspinning and 
a rescaling. In the present example, the rescaling is fixed by the requirement of keeping 
H unchanged in the infinite volume limit. More specifically, if we define 

^<max-l'-'^i ^ ]j^i(t^Xn^a.,-,^2,l + ,...,X2,o) (3-5) 

with 

^rimax-i ~ ^rimaxl ■ ■ ■ — X2 , (3.6) 

we can rewrite the energy as 
with 

^f' = -2 E (j)" , E , (E (3-8) 

In other words, we can blockspin without thinking about H and then include the second 
term of equation (3.7) in the local measures. The problem is then identical to the original 
problem except for the fact that we have to reduce rimax by 1 and to use a new local 

measure. 

The change in the local measure can be expressed through the recursion relation 

>V„+i(0') = 7V„+i exp((/3/2)0'2) /"^" de>V„(^ + 0>Vn(^ - , (3.9) 

where Nn+i is a normalization factor which can be fixed at our convenience. The symbol 
yy has been used instead of W in order to specify that the field had been rescaled at 
every step. The exphcit relationship is 

>V,((c/4)"/Vn) = W^n(0n) , (3.10) 
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where 0„ denotes the sum of all the fields in a block of size 2". Introducing the Fourier 
representation 

^e^'-^R^ik) , (3.11) 

the recursion formula becomes 

Rr,+,{k) = C„+i e-^^^ [RniV^ , (3.12) 

with Cn+i another arbitrary normalization constant related to the previous one by the 
relation C„+i = (c/-\/2)A^„+i. We will fix the normalization constant C„ is such a way 
that Rn{0) — 1- Rn{k) has then a direct probabilistic interpretation. If we call 0„ the 
total unrescaled field inside blocks of side 2" and < ... the average calculated 

without taking into account the interactions of level strictly larger than n, we can write 

RM = 1 + tl^^c/iy-/' < (0„)^ . (3.13) 

We see that the Fourier transform of the local measure after n iterations generates the 
zero-momentum Green's functions calculated with 2" sites. 

Everything that has been done in this section can be generalized in a straightforward 
manner for models with N components. All we need to do is to replace by a iV- 
dimensional vector (p, k by k and by d^^. The N component model is discussed in 
section 13. 

3.3. The Gaussian (UV) fixed point 

In the case where the initial measure is Gaussian 

= >Vo(0) = e-^-^' , (3.14) 

we have 

W„(0) oc e""^""^' (3.15) 
with Aq = A and the An calculable from the recursion relation 

A„+i = -/3/2 + (2/c)A„ (3.16) 
which follows from equation (3.9). The fixed point of this transformation is 

With our assumption < c < 2, A* is positive and finite. After resumming the terms, 
we obtain 

A„ = A* + (A-A*)(2/c)" . (3.18) 
Using equation (3.10), we then find that the measure for the main field $ is 

WnmaA^) oc e-^"— *'(^/^)"-"^ . (3.19) 
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This implies that 

K//(0c) = {{A - A*) + A\c/2Y-^)(t>l . (3.20) 

Given that < c < 2, the second term disappears in the infinite volume limit and if 
A> A* can interpret A — A* as a. quantity proportional to the square of the mass. 
On the other hand if A = ^4*, we have a massless theory in the infinite volume limit. 
At finite volume, we expect that for A — A*, 

oc (c/2)"™"- oc , (3.21) 

with L the linear size of the whole system defined by 

= O = 2"™"^ . (3.22) 

Putting these two equations together, we obtain 

c = 2^-2/-° . (3.23) 

In other words, the parameter c can be tuned in such a way that a Gaussian massless 
field scales with the number of sites in the same way as a D dimensional model. 

Following the same reasoning, we can determine the parameters introduced in 
subsection 2.2. Since we integrate over two sites in one RG transformation we have 

= 2 . (3.24) 

On the other hand, comparing the field rescaling, we obtain 

(4/c) = (3.25) 
which is consistent with equation (3.23) only if 

77 = . (3.26) 
In summary, we have 

c = £^-2 . (3.27) 

Wc can reinterpret equation (3.18) by noticing that it implies that at each RG 
transformation, An — A* is rescaled by a factor 2/c. If we introduce some reduced 
variable pn = {An — A*)/ A*, then 

Pn+l = {2/c)pn = ^Vn ■ (3.28) 

With the notations introduced in equation (2.28), we could say that the exponent y at 
the Gaussian fixed point associated with the mass term squared is 2. 

It should also be said that the Gaussian fixed point is often called the UV fixed 
point because, it is common to use perturbation theory about the Gaussian fixed point 
to construct a RG fiow from the Gaussian fixed point to the nontrivial fixed point (which 
is then called an IR fixed point) . 
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3.4. The HT fixed point 

In the high temperature (HT) hmit, /9 = 0, equation (3.12) becomes 

Rn+iik) = (i?„(v^ k)Y . (3.29) 
We set Cn+i — 1 and require -R(O) = 1. It is easy to check that 

is a fixed point of equation (3.29) for arbitrary B. However if S ^ and the exponent 
of k is not a positive integer, we have a branch cut at k = and by taking sufficiently 
many derivatives with respect to k, we obtain expressions which blow up at A; = 0. For 
the range of values < c < 2, the only way to get a power of k in the exponential 
which is an integer is to have c = 1, which according to equation (3.23) corresponds to 
D = 2. For c = 1, the exponent of k is 1, and the Fourier transform in equation (3.11) 
is ill defined unless we replace k by \k\ which in turn leads to singular derivatives at 
k — 0. Consequently, the only choice that leads to a probability distribution with finite 
moments is B — 0, or in other words, R — 1. This fixed point remains a fixed point 
when /3 7^ and is called the HT fixed point hereafter. The HT fixed point corresponds 
to an arbitrarily narrow probability distribution about for the main field and can be 
interpreted as the case of an arbitrarily massive free field. 



4. Equivalent forms of the recursion formulas 

The recursion formula will be the main tool used for calculations hereafter. It is 
important to identify equivalent or inequivalent forms and to find accurate numerical 
implementations. In this section, we review equivalent forms of the recursion formulas 
used by Baker [6], Felder [13] and Koch and Wittwer [49, 50]. This section follows closely 
reference [50] but with different notations. In the following, the constants A^^, A^^', . . . 
need to be fixed by some additional requirement and play no essential role. 



4.1. Baker 's form 

It is sometimes convenient to factor out the Gaussian fixed point so that in a new 
"system of coordinates" , the Gaussian fixed point is represented by a constant. If we 
define 

}Vn{4>) = exp(-^*02)7i„(0) , (4.1) 
the recursion formula becomes 

Hn+M = K^, / d^CXp{-2A*e) nn{^ + 0^n(^ ' , (4.2) 

J-oc Vc vc 

It is clear from this equation that if 7i„ is a constant then Tin+i is also constant. What is 
remarkable about this redefinition is that we have replaced exp((/9/2)0^)(i^ in equation 
(3.9) by d^exp{—2A*^'^). Any other multiplicative redefinition, would in general lead 
to some hybrid form where the weight depends on (p and ^. 
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It should also be noted that the variance of the Gaussian weight can be replaced 
by another value by a simple change of variable. If we define 

n^^K<l>)^n(B<t>) , (4.3) 

the recursion formula becomes 

J-oo 

The recursion used by Baker in reference [6] is obtained by setting 

^p/2A^)V^](^) = exp(-ig„(0)) (4.5) 

Note that Baker allowed a non-zero ri in order to mimic the scaling near the non- 
trivial fixed point instead of the scaling near the Gaussian fixed point for a conventional 
lattice model. It is possible to take into account this modification by changing the 
value of D used in equation (3.23). Confusion can be avoided by making reference to 
the rescaling factor c~^/^. More exphcitly, if we call D' the dimension used for a field 
rescaling that depends on rj, then the correspondence is 

2{2-n-D')/2D' ^ ^-1/2 ^ 2^2-D)/2D _ j^^ g^ 

A short calculation shows that D' = D{1 — {ti/2)). As an example, if we want to have 
a scaling corresponding io D' — ?> and 77 = 0.04, we can simply work with D — 3.061..., 
with D defined in equation (3.23). Using notation that we hope make clear the relation 
with reference [6] , the recursion formula can then be written as 

exp(-ig„+i(0)) = N';^, jy exp Ke - \Qn{^^ + - ^Qn(^ - O) (4.7) 
4. . 2. Gallavotti 's form 

An alternate way of formulating the recursion is to use a convolution with the Gaussian 
fixed point: 

/+00 
(i0cxp(-A'^(t-0)2)H„(0) . (4.8) 
■00 

By factoring out the two exponentials of the quadratic terms, it is possible to relate Q 
to the Fourier transform of W introduced in equation (3.12). More explicitly: 

ex-p{A*t'^)gn{t) oc Rn{i2tA*) , (4.9) 

and in particular the HT fixed point is now proportional to exp{—AH^) Under the new 
transformation, the recursion formula takes the form: 

^ d^0M-2j-^e) [Sni^^+O) . (4.10) 

This form of the recursion formula was used in references [49, 13, 50] and its origin 
can be found in the work of Gallavotti [51]. More about this question can be found in 
subsection 5.2 below. In the following, we call this form of the recursion formula the 
Gallavotti's form. 
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the recursion formula becomes after some algebra: 

From the definitions of Qn and Tm we have 

Tn{t) (X Rn{i2tA*) , 
and we can prove the equivalence of (4.12) and (3.12) by using the identity 
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Finally, if we write 



"+00 
X / exp 



fn{t)^Tn(t 



2-C 



A* {A - c) 

we obtain the form most often used in reference [50] : 

1 



/ 

J — c 



ds exp( 



1-c/A 



4-3. Summary 



(4.11) 



(4.12) 



(4.13) 



(4.14) 



(4.15) 



(/„(s + iv^)) . (4.16) 



The form of the U V and HT fixed points in the various set of coordinates is summarized 
in table 1. 

Table 1. Form of the fixed point up to an overall constant in the various coordinates. 
1 is short for a constant and S short for a delta function; the last line is the defining 
equation. 





W 


n 


G 




R 


/ 


uv 




1 


1 




C 4A* 


2-c .2 


HT 


6 


6 




1 


1 


1 


Eq. 


3.9 


4.1 


4.8 


4.11 


3.11 


4.15 



All this looks quite reminiscent of quantum mechanics where we can look at 
a problem in a basis where the position operator is diagonal or in another basis 
where the momentum operator is diagonal. The two basis are related by a unitary 
transformation. More generally, in quantum mechanics, unitary transformations do not 
affect the spectrum of hermitian operators (which represent observables) . Universality 
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is a stronger notion than observability in the sense that some non-universal quantities 
may have an absolute physical meaning independent of our choice of integration 
variables in the functional integral. The idea of transformations that leave universal 
properties unchanged have been discussed in reference [52] where this translate into a 
reparametrization invariance. 



5. Inequivalent extensions of the recursion formula 

In this section, we discuss more general recursion formulas which coincide with the HM's 
one for a particular choice of the parameter £ that controls the change in the linear scale 
after one RG transformation. 



5.1. Relation with Wilson's approximate recursion formula 

Wilson's approximate recursion formula was the first simplified RG transformation that 
was proposed. It appears in one of the basic RG papers [3] as a result of a rather 
involved analysis of the partition function of a scalar model with a UV cutoff. It played 
an important role, because it was realized that the RG method could lead to a fully 
numerical treatment without any reference to expansions such as perturbation in a weak 
coupling. It is intended to represent a situation where i — 2, but instead of having 2^ — 1 
integration variables, it has only one corresponding to an approximation where the 2^ 
fields take only two independent values [32], one in each half of the block. Having 
decoupled the number of integration variables to i, we can now write for arbitrary i: 

d^e-^' fel (^^--0 + d Hi:^ (e-^ -A] \ (5.1) 

oo 

For £ = 2, we obtain Wilson's approximate recursion formula, most often written in a 
form similar to equation (4.7). For i — 2^/^, we obtain the HM recursion formula in 
the form given in equation (4.2). 

The general case £ = 2'^ was discussed in reference [53] where it was shown that for 
D = 3, as ( increases from ^ = 1/3 to 1, the exponent 7 decreases monotonically from 
1.30 to 1.22. Clearly, different values of C correspond to different classes of universahty. 



5.2. Gallavotti's recursion formula 

The recursion formulas presented in subsection 4.2 can also be extended for arbitrary 
i. This is indeed easier because the number of sites integrated for the HM, namely 2, 
appears as the exponent for Q and /. The replacements are 

2 (5.2) 

7 ^r^-^ (5.4) 

4 
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For instance equation (4.16) becomes 



/•+00 



1 



s') [fl^\s+tr^/'-')f .{5.5) 



ds exp 



1 - ^-^-2 



Again for i = 2V-D 

we recover the HM recursion formula in the form given in equation 
(4.16). In reference [51], Gallavotti has introduced equations that can be identified with 
the case £ = 2. The variable s in equation (5.5) plays a role similar to the Gaussian 
variable in his notations, and the factor = 2^ is reabsorbed in the definition of 
the potential. The limit £ — > 1 of equation (4.10) will be discussed in section 9. 

For £ = 2 and D = 3, the value 7 = 1.30033... was obtained numerically in 
reference [54]. This value is significantly different from the value 7 = 1.2991407... 
obtained in references [55, 56] for £ — 2^/^ and D — 3. For D — 3, the limit £ — > 1 
was studied in references [57, 17] with the result 7 = 1.299124. The difference in the 
fifth digit is significant and was confirmed by new calculations [19] (see also our section 
9 below). This shows that again different values of £ correspond to different classes of 
universality and also that for a given i ^ 2^/^, the extensions given in equations (5.1) 
and (5.5) are inequivalent. In addition we see that the i dependence is much weaker 
in equation (5.5) and the slope is apparently opposite to the slope found for equation 
(5.1). We are lacking calculations at intermediate values but as far as we can see, when 
£ increases, 7 increases. 

6. Motivations, rigorous and numerical results 
6.1. Motivations 

Dyson's hierarchical model was invented and reinvented several times with different 
motivations that we briefiy review. Dyson's original motivation [5] was to construct 
models more weakly coupled than the one dimensional Ising models with long range 
hamiltonians of the form 



Dyson's was trying to figure out if the model has an ordered phase when a = 2. Dyson 
constructed a more general family of models where c' is replaced by bi in equation (3.1). 
He proved several theorems concerning the infinite volume limit and the existence of 
phase transitions at finite temperature that are discussed in the next subsection. 

The fact that o; = 2 is borderline can be anticipated by refining [58] Landau's 
argument for the absence of ordered phase in one dimension. If ck < 2 and a system of 
size L has an average magnetization /i > 0, then the cost in energy for hipping all the 
spins in a large subsystem, is of the order of L^"". If a > 2, the cost should grow slower 
than ln(L) when L is increased. On the other hand, there is of the order of L ways 
to choose the subsystem and so the gain in entropy is of order kTlnL. Consequently 
if q; > 2, the entropy dominates and the free energy is decreased after flipping, which 
is incompatible with the possibility of an equilibrium situation. On the other hand if 




(6.1) 
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a < 2, there is no incompatibility. In the case a = 2, both contributions are logarithmic 
and a more careful estimate is required. Thoulcss [58] estimated that if the width of 
the magnetization distribution is proportional to L^/^, the change in free energy is 

AE = 2//VlnL - (l/2)A;rinL , (6.2) 

and an ordered state seems possible for T small enough. If this occurs at some strictly 
positive critical temperature Tc, then the magnetization changes abruptly to {kTc/^JY^'^ 
when the temperature is lowered to Tc- This is called the Thouless effect. Later, the 
existence of an ordered state at sufficiently low temperature for the model defined by 
equation (6.1) with a = 2 was proved rigorously [59] as well as the Thouless effect [60]. 

Baker rediscovered the HM [6] in the context of the development of the RG ideas. 
His goal was to construct models for which a simple recursion formula would be exact. 
He reinvented Dyson's HM and several variant of it. This has been partially reviewed 
in sections 4 and 5. 

The hierarchical structure of the block variables can be naturally reconstructed 
using the 2-adic numbers (see section 14). At the end of the eighties, physicists started 
reformulating models of classical, quantum and statistical mechanics over the fields of 
p-adic numbers [61]. In particular models of random walks over the p-adic numbers 
were considered [62, 63, 64] and it was recognized that it was possible to reformulate 
the HM as a scalar model on the 2-adic fractions [65]. This reformulation was used 
in high temperature (HT) expansions [66], helps understanding the absence of certain 
Feynman graphs [3] and suggests ways to improve the hierarchical approximation [21]. 
This is discussed in more detail in section 14. 

6.2. Rigorous results 

Dyson proved several theorems for the HM with 

c = 2^-" , (6.3) 

where a is the same as in equation (6.1). With this notation it is clear that the 
ferromagnetic interactions arc weaker than for the model of equation (6.1). If we 
refer to figure 1, the relative strength of the couplings between the 0th spin and 
its right neighbors are, from left to right, 1, 2~", 3~", . . . , 7~" for equation (6.1) and 
2-a^4-a 4-a g-a g-a g-a g-a f^j. ^j^g fj|^ Griffiths has proved that for Ising models 

with ferromagnetic interactions, the averages of two arbitrary spins variables are positive 
and increase with the strength of the ferromagnetic interactions. Gonsequently, if 
one can prove that < 0„ Ylim > blows up in the infinite volume limit below some 
temperature for the more weakly coupled model, it will also blow up for the other model. 
Dyson proved the existence of the infinite volume limit for a > 1 and that there was a 
phase transition at finite temperature if and only if 1 < a < 2. In other words, there is 
no phase transition at a = 2 for the HM, but this does not allow us to extend the result 
to the more strongly coupled model of equation (6.1) that has indeed a phase transition 
at finite temperature [59]. With the notations used in section 3, Dyson theorems mean 
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that an infinite volume Umit exists for c < 2 (D > 0) and that a phase transition at 
finite temperature occurs if 1 < c < 2 {D > 2). 

The HM has been studied near using the e expansion. The existence 

of a non-trivial fixed point for e small enough was proved [67]. The e-expansion was 
shown to be asymptotic [68]. Many details regarding this approach can be found in 
reference [30] . Extensions beyond the hierarchical model are discussed in reference [69] . 

The HM has also been studied directly a,t D — 3 {c — 2^/^). The existence of a 
non-trivial fixed point was proved for a large enough number of components [70] . Proofs 
of the existence of the non-trivial fixed point in D = 3 for the one component model 
were given in references [49, 54, 50, 71, 72]. They also put exponential bounds on the 
fixed point in various representations discussed in section 4. In particular, for real 0, 
some positive C and 7i* the non trivial fixed point of equation (4.2), the following bound 
holds: 

|7Y*(0)| <exp(-C(/)^) 

6.3. Numerical results 

The critical and tricritical behavior of the HM was investigated numerically in the 
presence of a magnetic field and a staggered magnetic field [73] . The literature contains 
many numerical estimates of the critical exponent 7 for D = 3. The first calculation 
was done by Wilson. The result reported in reference [6] is 7=1.2991. In the following, 
unless specified difi^erently, errors of order one should be assumed for the last printed 
digit. Values when 2/D is a multiple of 0.05 are given in reference [74]. Interpolating 
linearly to D = 3, we obtained 1.302. Using the e expansion up to order 34 and a Borel 
resummation method, the value 1.2986 was obtained in reference [75]. Analysis of the 
HT expansion [27, 28] yields 1.300(2) for D ^ 3. The value 1.299141 can be obtained 
from a footnote in reference [54]. The value 1.29914 was obtained in reference [76]. 
Using two independent methods discussed in section 7, the value 1.299140730159 was 
found in references [55, 56]. Less accurate calculations in the low-temperature phase 
were performed in reference [77] and confirmed hyperscaling with three decimal points. 

7. Numerical implementation 
7.1. Polynomial truncations 

The recursion formula can be implemented numerically using numerical integration 
methods in equations (4.2) or (5.1) as was done for instance in Refs [3, 73, 53]. 
However, in the symmetric phase, it seems easier to get very accurate results by using 
polynomial approximations in forms based on the Fourier transform such as equations 
(3.12) or (4.16). This method was justified rigorously in [54, 72], and used for numerical 
calculations for instance in references [76, 78, 56, 14] . For definiteness, we start with the 



(6.4) 
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recursion formula for R given in equation (3.12). A finite dimensional approximations 
of degree l^ax has the form: 

Rn{k)^l + an,ik^ + an,2k^ + ... + a„,z_A;2'— . (7.1) 

This type of approximation can be justified in the context of the HT expansion and 
works extremely well in the symmetric phase. We can reabsorb the inverse temperature 
in /c in such a way that it does not appear anymore in the exponential in equation 
(3.12). This change would then be compensated by a transformation ^ f3''an,i and 
the truncation at order A;2''"°^ would be sufficient to calculate exactly the HT expansion 
of R up to order /J'™"^ . This technique was used in references [27, 28] . It was then realized 
(by accident) that large order coefficients in the HT expansion could be calculated in 
good approximation by using polynomial truncations at order 10 times smaller than the 
HT order. The same method can be applied to numerical calculations in the symmetric 
phase. The apparent convergence is studied empirically in reference [78]. 

With the polynomial truncation, the recursion formula equation (3.12) becomes a 
/„i(j^-dimensional quadratic map. After squaring R, we obtain a polynomial of order 
2lmax in k'^- We could in principle truncate at order Imax, however, the derivatives in the 
exponential in equation (3.12) will lower the degree and the terms of order larger than 
Imax will contribute to the orders smaller than Imax after enough derivatives are applied. 
Of course, a truncation at order l^ax is made after all the derivatives are performed, 
but it was realized empirically [79] that intermediate truncations reduce the accuracy 
of the calculation. The explicit algebraic transformation reads. 

_ E?^r(Ep+,=^«n,pan,.)[(20!/(/ -m)!(2m)!](c/4)'[-(l/2)/3]^- 

"""""''"^ E£r(Ep+,=.«n,pan,.)[(2/)!//!](c/4)'[-(l/2)/3]' 

The initial conditions for the Ising measure is Ro{k) = cos{k). For the LG measure, 
the coefficients in the /c-expansion need to be evaluated numerically. The susceptibility 
at finite volume and higher moments can then be obtained by rescaling the coefficients, 
for instance: 

X(f) = -2a„,i(2/c)- . (7.3) 
7.2. Volume effects in the symmetric phase 

When calculating the susceptibihty at values of (3 slightly below /3c, we spend about 
— ln(/3c — (5)/ln{\i) iterations near the fixed point. During these iterations, the round- 
off errors are amphfied along the unstable direction (see next subsection). After that, the 

order of magnitude of the susceptibility stabilizes and the corrections get smaller by a 
factor I at each iterations. At some point, all the recorded digits stabilize (irrespectively 
of the numerical errors which occurred in the first stage described above). This gives 
the estimate [78] for the number of iterations n(/3, P) to stabilize P digits (in decimal 
notations) 
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7.3. The Eigenvalues of the Linearized RG Transformation 

The critical exponents can be calculated by linearizing the RG transformation near the 
fixed point R*{k) specified by the coefficients a*i. We express the coefficients after n 
iterations in terms of small variations about the fixed point: 

an,i = a* + ^0-11,1 ■ (7.5) 
At the next iteration, we obtain the linear variations 

Imax 

San+1,1 = ^ Mi^^5an,m ■ (7.6) 

m=l 

The l^ax X ^max niatrix appearing in this equation is 

Ml,m = ^ , (7.7) 

evaluated at the fixed point. 

Approximate fixed points can be found by approaching f3c from below and iterating 
until the ratio a„+i,i/a„,i takes a value which is as close as possible to 1. The 
determination of Pc can be done by following the bifurcations in a„+i,i/a„,i for 
sufficiently large n. When /3 < /3c, the susceptibihty stabilizes at a finite value without 
subtraction and for n large enough, On+i.i/on,! c/2. On the other hand if /3 > /3c, the 
unsubtracted susceptibility grows like the volume and for n large enough, an+i,i/fln,i c 
(until the polynomial truncation breaks down). 

The approximated fixed points obtained with this procedure depend on P^- Using 
their explicit form which we denote R*{k, jSc), one obtains a universal function 9^(A;) by 
absorbing (3 into k 

^{k)^R\^/p,k,(5,) . (7.8) 

It was shown that in very good approximation, 9^(A;) is independent of the initial measure 
considered [55, 56]. Numerically, 

^{k) = 1. - 0.35871 134988A;^ + 0.0535372882A;^ - . . . . (7.9) 

This function is related to the fixed point fRwi^s^) constructed in reference [72] as 
follows by the relation 

m) « fKw{{^)k^) . (7.10) 
zc 

Extremely accurate values of the Taylor coefficients of fxw can be found in the file 
approx.t in [72]. The constant of proportionality is fixed by the condition 9^(0) = 1. 
The relation with the non-trivial fixed point /* of equation (4.16) is f*{it) = fKwi—t"^)- 

The first six eigenvalues of M/^m, from reference [56] are given in table 2. 
Using 7 = ln(2/c)/lnAi and = — lnA2/lnAi from section 2, we obtain 

7 = 1.299140730159 (7.11) 
A, = 0.425946858988 . (7.12) 
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// 


A. 


1 


1.42717247817759 


2 


0.859411649182006 


3 


0.479637305387532 


4 


0.255127961414034 


5 


0.131035246260843 


6 


0.0654884931298533 



These estimates of the exponents were in agreement [56] with those obtained from fits 
of numerical data near criticality based on the parametrization of equation (2.26) with 
12 significant digits for 7 and 6 significant digits for A^. These fits also provide the 
non- universal amplitudes. 

7.4- The critical potential 

It is possible to Fourier transform numerically '^{k) for a not too large value of the 
conjugate variable (p. The result is that apparently, the Fourier transform is a positive 
bell-shaped function with only one maximum at (f) — 0. We define 

= -ln^(0) , (7.13) 

and we obtain an apparently convex function with only one minimum shown as the solid 
line in figure 2. In order to compare with calculations done with the 7i„ formulation, 
we need to subtract A*(f)'^ from this function. Since we have removed the f3 dependence 
by reabsorbing it in the definition of k, we must calculate A* with (3=1. Subtracting 
(c/(2(2 — c)))0^ ~ 0.8510^ from the rescaled potential then becomes a double-well 
potential with minima at ~ ±1.688. This is illustrated in figure 2. This explains that 
figures with a minimum away from the origin appear in references [32, 73]. 

It should be emphasized that U((f)) is a rescaled potential and not the effective 
potential. The relationship between the two can be worked out from equations (3.10) 
and (3.11). Assuming we are exactly at the non trivial fixed point, W remains 9^(0). 
One then obtains 

Z oc i?*(-iJ(4/c)"""'-/2) , (7.14) 

and 

QK//(0c) = [/(0eC"-"-/') • (7.15) 

This result can be interpreted in the following way: at the fixed point, the only scale 
in the problem is the size of the system and all the quantities scale with it according 
to their dimension. If we perform Taylor expansions in equations (7.14) and (7.15), we 
see that when rimax the coefficients of Z blow up and the coefficients of (fr^^ 
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blow up if r > D/{D — 2). In the same limit, the coefficients of 0^'' tend to zero if 
r < D/{D — 2). The limiting case corresponds to Vc = D / {D — 2) where the coefficient 
is volume independent. For instance for D = 4, Vc = 2 and for D = 3, Tc = 3 which 
corresponds in both cases to the marginal direction in perturbation theory. It should 
however be noted that if we work at finite 4>c, we are in fact probing U at large value of 
its argument and the Taylor expansion may not be convergent (the radius of convergence 
should be determined by the complex zero of 9^(0) closest to the origin). The above 
determination of Tc implies that, 

C/(0) oc = l^l^/^" for|(/.| ^ oo . (7.16) 

For D = 3, the (f)^ behavior is compatible with the rigorous bound of equation (6.4). 

Equation (7.16) can be obtained from the scaling hypothesis discussed at the end 
of section 2.2. First, we find from equation (7.14) that 

i^G[J] = - ln(i?*(-iJ^(^+2)/2^^ ^ constant , (7.17) 
which implies that for large values of | J| and \4>\ 

G[J] oc I J|2^/(^+2) ^ I j|D/y, ^ (7_;^g) 

and 

101 oc |J|(^-2)/(D+2) _ (7^^g) 

Combining the two above equations and the Legendre transform, we recover equation 
(7.16). 

7.5. The low temperature phase 

The calculation in the low temperature phase requires the introduction of a constant 
external magnetic field coupling linearly to Since 

W„{<Pn,H)^W„{<Pn)e'''^'' , (7.20) 
after Fourier transforming and rescaling one obtains [77] 

The connected Green's functions can be obtained by taking the logarithm of this 
generating function. The average are understood at non-zero H. In order to observe 
the magnetization, it is essential to take the infinite volume limit before taking the limit 
— > 0. For any non-zero if, no matter how small its absolute value is, one can always 
find a n large enough to have \H{A/cYl'^\ > 1. The non-linear effects are then important 
and linearization does not apply. It was checked [77] that when such a n is reached, the 
value of the x^'^^ stabilizes at an exponential rate. One can then, first extrapolate at 
infinite volume for a given magnetic field, and then reduce the magnetic field in order to 
extrapolate a sequence of infinite volume limits with decreasing magnetic field, toward 
zero magnetic field. 
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In the following, we use the notation 

x(^)oc(/3-/3e)-^'' , (7.22) 

for the leading exponent. It is customary to use the notation 71 = — /3, but we avoided 
it here because of a possible confusion with the inverse temperature. Obviously 72 = 7- 
In the symmetric phase and for q even, we have the order of magnitude estimate 

^(9) ^ 2-"*(4/c)''"*/2 . (7.23) 

with n* defined as in section 2 by the relation — /^dA" = 1. Eliminating n* and using 
the expression of 7 as in section 2, we obtain 

7, = 7[(?/2)ln(4/c) - ln2]/ln(2/c) . (7.24) 

In the case D — 3 {c — 2^/^), this becomes 

-fg = 1.29914073 . . . X (5g - 6)/4 . (7.25) 

We will show below that this equation is a particular case of a more general relation that 
follows from the scaling hypothesis. In reference [77], the following numerical results 
were obtained 

7i= - 0.3247 

72 = 1.2997 (7.26) 

73 = 2.9237 , 



Nonlinear Aspects of the RG Flows of Dyson's Hierarchical Model 



32 



which agrees with three significant digits with the prediction of equation (7.25). 

Equation (7.24) follows from a slightly stronger form of the scaling hypothesis. The 
basic scaling relation of equation (2.27) is satisfied if we further assume that 

G,{t,J)^t^''g{J/t^'') , (7.27) 

for a well-behaved function g and with the gap exponent 

A, . I = (£+2-^ p^^, 

By construction the argument of g is invariant under the rescaling of equation (2.27). 
Each derivative with respect to J brings down a factor t~^'>. This implies that 

7, = -Du + qAg . (7.29) 

In the case of the HM, u — 7/2 and Ag = {D + 2)7/4 and we recover equation (7.25) 
for D — 3. For the HM in arbitrary dimension, we could also write 

T, ^ -^^Sm^ . (,30) 

7. 6. Practical aspects of the hierarchy problem 

In absence of wave function renormalization, the square of the renormalized mass ruR 
in units of the UV cutoff A can be defined as the inverse susceptibility. Keeping the 
mass small when the cutoff increases requires a large susceptibility. In the calculations 
discussed above, a large susceptibility is obtained by fine-tuning /3. However, we can 
also keep (3=1 and fine tune another parameter such as the bare mass me in a Landau- 
Ginzburg potential. In this case, we have 

mVA ~ {\mBc - mB\/KV'^ , (7.31) 

In four dimensions, 7 = 1 and if we take m = 100 GeV, a typical electroweak scale, and 
A = 10^^ GeV of the order of the Planck mass, we need to fine-tune with 34 digits. 
This is often called the hierarchy problem and seen as an argument against fundamental 
scalars [80]. The main virtue of the RG approach is to separate the relevant and 
irrelevant part of the information contained in the partition function. At each iteration, 
the information relevant to understand the large distance behavior is amplified, while 
the rest of the information is discarded according to its degree of irrelevance. However 
if some "noise" is introduced in this process, for instance as round-off errors in the 
calculation, the error in the relevant direction will be amplified too. This may lead to 
situations where the amplified errors wipe out the final result. In the case of the HM, 
the problem can be solved by increasing the arithmetic precision in the implementation 
of equation (7.2). This is documented in reference [81]. 
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8. Perturbation theory with a large field cutoff^ 

8.1. Feynman rules and numerical perturbation theory 

An attractive feature of the HM is that it is possible to calculate perturbative series to 
large order by blockspinning numerically, order by order in an expansion parameter, for 
instance A for a perturbation. This method can be used analytically to reconstruct 
the Feynman rules [82, 16]. In practice, the diagrammatic expansion is much more 
complicated than the numerical method. However, for comparison with calculations 
based on diagrams, it is useful to know the Feynman rules. 
For an initial measure of the form 

Wo{<P) = exp ( - {A^ + im|)02 _ A^^) , (8.1) 
we obtain the usual Feynman rules for a theory with the following replacements: 



{2n) 



D 



J] 2-"-^ (8.2) 



n=0 

1 1 



with k'^{n) = 2A*(c/2)". The interpretation is quite simple, the integral over the 
momenta is replaced by a sum over momentum shells similar to those introduced by 
Wilson [3, 32]. After one RG transformation the UV cutoff A is lowered to A/i = 2"^/-^A 
and the volume of momentum space in D dimension is reduced by a factor 2. The volume 
of the 0-th sheU is 1/2, the volume of the 1st sheU is 1/4 etc... Similarly, (c/2)" = ^-^n 
represents the square of the momentum in the n-th shell. 



8.2. Perturbation theory with a large field cutoff 

It is well known [83] that perturbative series are in general divergent. Their zero radius 
of convergence is due to large field configurations [84, 85]. However, the large field 
configurations have very little contributions to observables involving a few fields such 
as the magnetic susceptibility or the 4-point function. 

This point was realized by the author in two different circumstances. The first is 
quantum mechanics, quantum field theory in 1+0 dimensions, where the field variable 
is usually denoted x. A large field cutoff can be implemented by imposing that the 
potential becomes +oo at x = ±Xmax- If the field cutoff Xmax is large enough, the effects 
on the low energy levels are exponentially small [86, 87, 88]. The second circumstance 
is the HM [78] . The numerical procedure described in section 7 is based on polynomial 
approximations and is purely algebraic, however, we need to input Ro{k). So we need to 
do one integral numerically to start, namely the inverse Fourier transform of equation 
(3.11). At the end, we need to Fourier transform if we want to extract the effective 
potential. In doing the initial integral numerically it is convenient to introduce a large 
field cutoff and then monitor the effect of this cutoff when it is increased. It is clear 
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Figure 3. First nine perturbative coefficients (left to right) for the two-point function 
in unit of their infinite field cutoff value, as a function of the field cutoff (pmax- 



that for local measures that decay sufficiently fast, the effect is exponentially small for 
observables involving a few fields. 

On the other hand, the large order of the perturbative series involve averages of 
large powers of the field and is sensitive to the field cutoff. This is illustrated in figure 3 
where the perturbative coefficients of the zero-momentum two point function for = 3 
are plotted in units of their value at infinite field cutoff as a function of the field cutoff. 
One can see that for a fixed large field cutoff (pmax, some low order coefficients may be 
close to their asymptotic values, a few coefficient may be in the crossover region and most 
coefficients are much smaller by several order of magnitude than their asymptotic value. 
These three regimes are reminiscent of the three regimes encountered when calculating 
renormalization group flows between two fixed points [89, 90, 79, 38]. Note also that 
the shape of the transition seems universal as in the anharmonic oscillator case [88, 91]. 

This suggests a connection between the crossover observed in the behavior of the 
perturbative coefficients and the crossover behavior of the RG fiows. When we construct 
the RG fiows starting near the Gaussian fixed point and let them evolve toward the high- 
temperature fixed point, it should be possible to describe the first iterations using the 
Gaussian scahng variables (see section 10). On the other hand, after a large number 
of iterations, the scaling variables of the HT fixed point are the relevant ones. If we 
use regular perturbation theory, we expect that it will be impossible to find a region 
where the two expansions are valid due to the zero radius of convergence of the weak 
coupling expansions. On the other hand, if a field cutoff is introduced, the weak series 
have a nonzero radius of convergence and the direct calculations of critical amplitude 
as in reference [79] might be possible. A generic feature that we then expect is that 
if we calculate the perturbative coefficients with a field cutoff, by blockspinning, the 
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Figure 4. Number of significant digits obtained with regular perturbation theory at 
order 1, 3, 5, 15 (sohd, turning clockwise with order) and with (pmax = 2 (dash 
line), at order 1, 3, 11 (moving up with order) as a function of A, for the two-point 
function of the HM. 

first coefficients sliould stabilize quickly, while the large order in perturbation should 
stabilize after more iterations. This property was verified in reference [91]. 

8.3. Improved perturbative methods 

The field cutoff significantly alter the accuracy of the perturbative series. This is 
illustrated in figure 4 where the accuracy of perturbation theory for the two point 
function at various order is shown in regular perturbation theory and with a particular 
field cutoff. The figure makes clear that at sufficiently large coupling, the modified 
series becomes more accurate than the regular series. It is also clear that for a given 
field cutoff, the accuracy peaks near a specific region of the coupling. It is likely that 
at a given coupling, it is possible to find an optimal field cutoff that can be determined 
approximately using a strong coupling expansion as in a simple integral discussed in 
reference [92]. 

8.4- Large field cutoff in ERGE 

Understanding and controlling the large field configurations is an issue that goes beyond 
the scope of perturbation theory. In particular, it appears in the context of the RG flows 
of the effective actions [93]. It was noticed [94] that the introduction of a background 
field suppresses large field contributions to the fiows. We should also mention functional 
generalizations of the Callan-Symanzik equation [95, 96] where a running bare mass 
controls large fiuctuations. 
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9. Relation with the ERGE in the LPA approximation 

9.1. Polchinski equation in the LPA 

As explained in the first sections, the RG transformation of the HM can be reduced to a 
a simple integral equation because of the very special form of the non-local interactions. 
In general, the real space RG seems difficult and one may prefer a formulation in 
terms of the Fourier transform of the fields. An UV cutoff can be introduced and 
lowered in variety of ways (sharp cutoff or smooth cutoff functions, ...). An effective 
action can be obtained by lowering the cutoff, or varying a parameter in the cutoff 
function in such a way that the large momentum components of the fields get more 
integrated. The derivative of the effective action with respect to the cutoff (or a related 
parameter) can then be expressed in terms of an integral over the momenta of a function 
of the action and its derivatives with respect to the Fourier transform of the fields. 
This equation is exact and is often called an Exact Rcnormalization Group Equation 
(ERGE). This idea was introduced and developed in references [4, 97, 98, 11] and has 
generated a large interest that is still ongoing. Progress have been reviewed for instance 
in references [99, 100, 8, 9, 10]. An ERGE can be rewritten as an infinite set of coupled 
partial differential equations. A simple starting point is to neglect the evolution of the 
terms in the effective action involving derivatives. This is called the Local Potential 
Approximation (LPA). 

A simple equation that can be written in this approximation for an ERGE with a 
smooth cutoff function is [98]: 

- = W+(l--)0--(-)^ + ^ (9.1) 

where V(t, 0) is the effective potential and t a parameter that increases when the cutoff 
decreases. 

9.2. Infinitesimal form of Gallavotti's recursion formula 

We now consider the extension of Gallavotti's recursion formula equation (4.10). We 
introduce the notations: 

^„(0) = e-^(^'<^) (9.2) 
^„+i(0) = e-^(^/^"^) (9.3) 

In order to take the limit ^ — > 1, we need to fix the variance of the Gaussian weight and 
the overall normalization. We rewrite the extension of equation (4.10) for an arbitrary 
£ with an Gaussian weight parametrized in terms of a function p{£) as 

exp (-V(A/<.^)) = £°°^ exp (-|^-^«V(A. £'-/V+o) .(9.4) 

We require that p(l) = and p'(l) 7^ 0, otherwise p is arbitrary. In reference [13], we 
havep(£) =2(£-l). 
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The factor (7rp(^))~^/^ guarantees that when i ^ I, the two sides of equation (9.4) 
are equaL We now write £ ~ 1 + 5 and expand in 5. As mentioned, the terms of order 
cancel. Terms of order 5^/^ appear but become zero after integration over ^. Equating 
the terms of order 6, we obtain 

_A_ = W + (1 - -)«,- - [i- f - ^] . (9.5) 

If we write A = e^* and change the (p scale in order to get rid of the factor ^^-j^, we 
recover equation (9.1). 

It should be pointed out that despite the fact that the differential equation comes 
as the coefficient of the order S in the expansion of equation (9.4), equation (9.5) is 
not a linearization of equation (9.4). Indeed, when we reabsorb p in the integration 
variable ^, we introduce a term of order ^/p in the argument of the exponential. The 
terms that survive integration are quadratic in y/p, namely the second derivative of u 
and the obviously nonlinear square of the first derivative of u. This point was not fully 
understood in reference [18]. 

Note also that the Gaussian fixed point corresponds to = 0. The problem 
of finding the eigenvalues of the linearized RG transformation reduces to a time 
independent Schrodinger problem. The (pd/dip term can be eliminated by minimal 
substitution, introducing a (j? terms and the problem can be mapped into the problem 
of finding the eigenvalues of an harmonic oscillator. 



9.3. The critical exponents of Polchinski's equation 

A related equation is the so-called Polchinski equation. It can be written for N 
components as [101]. 

^ = -§u" + (l + - + {2-d)y- 2yu)u' + (2 - u)u , (9.6) 

— * — * 

with y — (f) . (f), u — 2V' and the prime denotes derivatives with respect to y. This 
equation can be derived [102, 101] from an ERGE due to Polchinski [11] using the LPA. 

For N=l, one can see that it follows from equation (9.1), by reexpressing it in terms 
of y and its derivatives and taking the derivative with respect to y of the resulting 
equation. Equation (9.1) is also obtained as the LPA of an ERGE due to Wilson [8]. 
The exponents were calculated in reference [101]. In particular for A^ = 1, they found 
7 = 1.2992 which is close to the HM value. A more precise value 7 = 1.29912 was 
obtained in reference [17]. 

In addition, Litim [15, 103, 57] proposed an optimized ERGE and suggested 
[104, 93] that it was equivalent to the Polchinski equation in the local potential 
approximation. The equivalence was subsequently proved by Morris [105]. The value 
of 7 for the optimized ERGE [57] in the case A^ = 1 is 1.299124 and differ by 2 in 
the 5th decimal from the HM. More recently [19], the calculations using the optimized 
ERGE and Polchinski equation were both repeated with more accuracy and compared. 
The numerical difference between the exponents of the two (analytically equivalent) 
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formulations was reduced to 10~^^. Their final result is 7 = 1.2991235477613 which 
confirms the non-equivalence with the HM. This question is also discussed for N > 1 in 
section 13. 

9.4- Infinitesimal form of Wilson approximate recursion formula 

The infinitesimal form of Wilson approximate recursion formula can be derived by 
following the same steps as for Gallavotti's recursion formula. First, we write Hn[(/)] = 
cxp{—Q{(f))). We then use the arbitrariness of the scale of the fiuctuations ^ as 
previously. The only difference is that the term of order a/5 disappear from 

Q((l + 5)1-^/^ + ^/SO + QUI + 5)'-^/^ - ^0 ■ (9-7) 

Consequently, there seems to be no {dQ/d(f))^ term in the final equation. This point 
was also noticed in reference [36]. This suggests that the limit £ — > 1 is Gaussian 
(7=1). In reference [53], a numerical calculation of 7 was done for values of £ = 2''. 
C = 1 corresponds to Wilson's case while ( = 1/3 corresponds to the HM. The limit 
we are interested in for the infinitesimal form is C ~^ 0. Unfortunately, in this limit, 
the numerical procedure used in reference [53] becomes unstable because the errors in 
the integration routine become more important as we need to iterate more times the 
basic recursion formula. Figure 1 in reference [53], indicates that 7 keeps increasing as 
C decreases. The last data point is for ( — 0.3. For smaller values of ( large errors bars 
develop as can be seen by repeating the calculation at closely chosen values of (. For 
instance near ( = 0.15, we found values of 7 as high as 1.34 and as low as 1.28. This 
calculation should be repeated with more accurate integration methods. 

9.5. Finite time singularities 

It has been argued that some ERGE in the LPA have finite time singularities [106] . This 
is not surprising given that the solutions of the fixed point equation for equation (9.6) are 
generically singular. More precisely, if we assume that du/dt = 0, we obtain a second- 
order differential equation for u. The solutions blow up at finite y for generic "initial" 
conditions at y = 0. This means that the derivative of the potential becomes singular at 
a finite value of the field. It has been shown numerically [101], that regular solutions can 
be obtained for special initial values using the shooting method and that these solutions 
correspond to the non-trivial fixed point obtained with other formulations. Rigorous 
results concerning the existence of global stationary solutions of equation (9.4) can be 
found in reference [13]. 

On the other hand, the possibility of having finite time singularities in Monte Carlo 
RG calculations for nearest neighbor models is controversial [107, 108]. More generally, 
these reviews question the existence of renormalized or effective Gibbs measures defined 
by certain RG procedures. 

For the HM with the polynomial approximation discussed in section 7, singularities 
after a finite number of iterations cannot appear (as they cannot appear for a finite 
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dimensional quadratic map). More generally, it seems possible to prove the boundedness 
of Tin or Qn (defined in section 4), for finite n, for a large class of initial functions. 

9. 6. Improvement of the LPA 

The improvement of the LPA for ERGE as a derivative expansion is a well-developed 
subject [100, 8, 9, 10]. However progress is still needed in order to get estimates of 
the exponents which can compete in accuracy with the best methods available [35]. 
It is possible that the basic differential equations for the effective potential and the 
coefficients of terms involving derivatives of the fields could be worked backward, at 
finite £, in order to produce a set of manageable coupled integral equations. The 
improvement of the hierarchical approximation is discussed with completely different 
methods in section 14. This could lead to £ — > 1 equations that could be in turn 
compared with the existing ones. We hope that some communication between the two 
approaches will be developed in the future. 

10. The nonlinear scaling fields 

10.1. General ideas and definitions 

In the study of ordinary differential equations, a standard method [109] to go beyond 
the linearized approximation near a fixed point, consists in constructing new coordinates 
where the equations become linear. In the context of the RG method, these new 
coordinates are called the nonlinear scaling variables (or scahng fields) and were first 
introduced by Wegner [31, 110]. 

10.2. The small denominator problem 

Rectification procedures are usually plagued with the "small denominator problem" 
initially encountered by Poincare in his study of perturbed integrable Hamiltonians [111]. 
In the RG case, this question needs to be discussed for each fixed point separately. To the 
best of our understanding, for the HM in D = 3, the problem can be completely avoided 
(but in a non obvious way) for the HT fixed point, it is not present for the nontrivial 
fixed point and it is essential to generate logarithmic corrections to the scaling laws near 
the Gaussian fixed point [31, 112]. 

In the rest of this section, wc will discuss in detail the case of the scaling variables 
associated with the HT fixed point and the nontrivial fixed point. Later we show that 
they can be combined in order to calculate non-universal critical amplitudes. The way 
the small denominator problem can be avoided for the scaling variables of the HT fixed 
point is interesting. At first sight, the construction seems impossible for D = 3 and more 
generally for rational values of because some of the denominators are exactly zero. 
A numerical study in D = 3 showed [113] that for all zero denominators considered, a 
zero numerator miraculously appears. Explicit calculations in arbitrary dimensions and 
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general arguments explaining why it should work to all orders were given in reference 
[114] and are summarized in subsection 10.5. 

10.3. The linear scaling variables of the HT fixed point 

As explained in section 7, the RG transformation in the symmetric phase can be 
approximated very accurately in terms of a quadratic map in a Imax dimensional space 

an+i,i = — , (10.1) 

with 

Un,a = rf a„,^a„,, , (10.2) 

and 

p.. ^ (-l/2)-^-'(2(^ + .))! 

" * ' ' (a< + !'-<t)!(2<t)! ' ^ ' 

ioY fi + u > a and zero otherwise. As in "relativistic" notations, the greek indices /i 
and u go from to Imax, while latin indices i, j go from 1 to Imax- Repeated indices 
mean summation unless specified differently. With the normalization of equation (10.1), 
a„, = 1 for any n and is not a dynamical variable. For small departure from the HT 
fixed point 5a„ ^ the linear RG transformation reads 

6an+i,i Mldanj , (10.4) 

with 

A.| = 2rf = 2(£rt-ir'^^J|^, (10.5) 

for i < j and zero otherwise. 

Ai is of upper triangular form and the spectrum is given by the diagonal elements: 

A(,) = 2{c/4y = £^--(^+2) , (10.6) 

in agreement with reference [30]. Note that the quantity in the exponent of i also appears 
in the exponent 7^ when g = 2r in equation (7.30). The tilde refers to the HT variables, 
eigenvalues etc... in order to avoid confusion with the same quantities for the nontrivial 
fixed point. As we assume c < 2 in order to have a well defined infinite volume limit (see 
section 6), all the eigenvalues are less than 1 and the fixed point is completely attractive. 
As r increases, the eigenvalues decrease and become more irrelevant. Equation (10.6) 
has a simple interpretation: the 2 stands for the volume increase and the (c/4)'' for the 
rescaling of the 2r— th power of the sum of the fields, which has a dimension a^/^"^^ in 
lattice spacing units, when the volume element is properly included, and of course 
assuming rj = 0. The eigenvalues in equation (10.6) can also be seen in equation (7.23) 
We call TZ the matrix of right eigenvectors: 

MX = A(,)7^[ , (10.7) 
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(with no summation over r). For convenience, the columns of TZ are ordered as the 
eigenvalues. We introduce the linear coordinates hi 

an,i = Tl\hn,r , (10.8) 
and which transform as 

hn+l,r — ^{r)hn,r (10.9) 

in the linear approximation. The matrix TZ\ and its inverse are also upper triangular and 
hi is of order j3\ just as a,„ / is. Wc fix the normalization of the right eigenvectors in TZ in 
such way that all the diagonal elements arc 1. This guarantees that hi = an,i + 0{l3''^^). 
In reference [114] it was proved that for the upper diagonal elements (j > i), 

and that 

(7^-l)i = i-iy-'T^i . (10.11) 

Note that because of the HT selection rules (i. e., the upper triangular form of the 
matrices), the matrix elements do not depend on the choice of Imax- 

It is easy to rewrite the exact RG transformation in the hi coordinates. Starting 
with the basic equation (10.1), we replace ao by 1 and ai by TZ^hp, we obtain a recursion 
formula of the form: 

~ ^ ~X^i)hi + A^%h^ 

1 + 2Al% + Al%h, ' ^ ' ' 

with coefficients calculable from equation (10.3). For instance. 

Pursuing the relativistic analogy, upper roman indices transform with TZ and the lower 
ones with {TZ)~^. 

10.4- The nonlinear scaling variables of the HT fixed point 

We now explain how to reexpress the linear variables hi in terms of the nonlinear scaling 
variables yi for which the approximate multiplicative transformation of equation (10.9) 
is assumed to be exact: 

yn+l,r = \r)yn,r ■ (10.13) 

If we use \n{yi) as our new coordinates, the RG flows become parallel straight lines. All 
the dynamics is then contained in the mapping. The monotonicity of these functions 
suggests a possible connection with field theory entropy [115], however the regularity 
near other fixed points may be an issue. 

As in reference [31, 110], we introduce the expansion 

hi^yi+ ^iMi..Sif2--- , (10.14) 

11,12, ••• 
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where the sums over the i's run from to infinity in each variable with at least two 
non-zero indices. In the following, wc use the notation i for (ii,i2, • • •)■ More generally, 
such vectors will be represented by boldface characters. The unknown coefficients Si^i 
in equation (10.14) are obtained by matching two expressions of /i+i,/, one obtained 
from the RG transformation of the hi given in equation (10.12), the other obtained 
by evolving the scaling variables according to the exact multiplicative transformation 
equation (10.13). The matching conditions can be written as: 

hn+l,l{K{y)) = hi{Xiyn,l, A2^n,2, • • •) • (10.15) 

and yield the conditions 

= ^ ■ (10.16) 



with 



j+k=i 



and 



+ E ^'J^^(m)AS^«P,k%r. (10.17) 
j+k+r=i TO 

A,i = A(,)-nA;-). (10.18) 



For a given set of indices i, we introduce the notation 

= Ytmm'^ . (10.19) 

m 

One sees that 2o is the degree of the associated product of scahng variables and Xi its 
order in the HT expansion (since yi is also of order /3'). Given that all the indices are 
positive and that at least one index is not zero, one can see that if j + k = i then 
2g(j) < Xg(i) and Xg(k) < Xq(i). Consequently, equation (10.17) yields a solution order 
by order in Xq or in Xi (since the r.h.s is always contains si^i of lower order in Iq and Xi) 
provided that none of the denominators Di ^ are exactly zero. 
We can now rewrite the denominators as 

Di,i = 2(^)' - 2^o«(^)^iW . (10.20) 

The parametrization of c in equation (3.23), implies that a a zero denominator appears 
when 

D-l{D + 2)^ DXo(i) -{D + 2)Xi(i) . (10.21) 

Given that the Xg are integers, this can only occur at some rational values of D. Ignoring 
temporarily this set of values, we can say that for generic values of c, the denominators 
are not zero. In the spirit of dimensional regularization, we can perform, order by order 
in the HT order Xi, the construction of the for a generic value of c and discuss the 
limit where c takes some special value at the end of the calculation. Since the linear 
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problem is completely solved and we may assume 2o(i) > 1. In addition, since both 
hi and yi arc of order f3\ we need Xi(i) > /. At lowest non-trivial order in we have 
X/(i) = /, and it has been shown [114] that, at that order, s;^i has only an apparent pole 
of order / at c = exactly canceled by a zero of the same order in the numerator. If 
Xi(i) > /, we can write 

Di,,^2{^y{c,rit.y-^^^'>Ti,,, (10.22) 



with 



and 



Ti,i = (c5«-' - c^^«-') , (10.23) 



Crit. = 4 X 2(i-^o«)/(^i»-') . (10.24) 

The only poles that we need to worry about are those where < Ccru < 2. An 
inspection [114] of the 175 terms up to order P'^ shows that all the poles at < c < 2 
were exactly canceled by zeros of the same order. Note that the maximally simplified 
rational expression for the do have poles but at values of c outside of the range 
< c < 2. A possible strategy for a proof would be to show inductively that for each 
term with < Ccrit < 2, the cancellation occurs and so the undesired poles do not 
propagate to higher order. Such an algebraic proof seems difficult because the i are 
rational expressions and the zeros at the numerator only appear after factorization of 
sums of such terms. It seems nevertheless reasonable to conjecture that si^i have no 
poles < Ccrit < 2. If this conjecture is correct, dimensional regularization provides a 
unique continuous expression for the coefficients for any c with < c < 2 and the model 
is formally "solvable" using the recursion for the coefficients given by equation (10.17). 
The conjecture implies that for any value of c in this interval, we can construct analytical 
expression of a„^; (which contains all the thermodynamical quantities) in terms of Oq,; 
(which depends on the initial energy density): 

an,i = {n-XhrCKVii^o), A^y2(ao), . . .) ■ (10.25) 

We will see that the initial values of y(ao) have a simple interpretation given in equation 
(10.41). 

It is also possible to express the nonlinear scaling variables in terms of the linear 
variables. Writing, 

yi-hi + Ylri,il[ht^ , (10.26) 

i m 

we can determine order by order the unknown coefficients r;,i of the expansion for generic 
values of c. 

10.5. Argument for the cancellation to all orders 

The generating function of the connected parts of the average values of the total field 
reads 

ln{Rn{k)) = al^^e + a^,2^^ + • • • , (10-27) 
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with 



a": 



J2 i-if°^^-\u^-iy-llfi- (10-28) 



i:Xi(i)=; m 



We are working in the HT phase and that we do not need to subtract powers of the 
magnetization. After a suitable rescahng of k described in subsection 7.1, we have 

= i-py^ilr^rn'r , (10.29) 

We assume that the initial values ao,/ are such that, 

hm xl?^ = x^^^ , (10.30) 

n— >oo 

is finite. In other words, we assume that P < Pc and that all the x^^^ are finite. As 
explained in section 6, this statement can be proved rigorously for a Ising measure. 
Prom equation (10.29), it is then clear that for n large enough, we have the leading 
scahng 

<zOc(2(^)')" = A^,). (10.31) 

This suggests a simple relationship between ; and ijn^i. Using Mobius inversion formula 
[116, 117], it has been shown that 

a^ = yi + 0{/3'^') . (10.32) 

Equation (10.32) means that there are no nonlinear contributions of order f3^ to a'(. For 
instance, there are no y\ or yiy2 terms in ag. This is expected because the nonlinear 
terms of order scale faster than y^, (assuming < c < 2). We say that a term "scale 
faster" , we mean that it goes to zero at a slower rate when n becomes large. In general, 
at each RG step, a term y'^ of order is multiphed by 

2^°(^)(|)'>A(o = 2(^)^ 

The strict inequahty comes from the fact that for the nonlinear terms 2o(i) > 1. It 
is thus clear that nonlinear terms of order f3^ would spoil the HT scaling of equation 
(10.31) and contradict the existence of a infinite volume limit. 

For higher order terms, the sign of the denominator Dn introduced in equation 
(10.18) tells us whether or not the term scales faster or slower than the linear term. 
With our sign convention, c > Ccrit.(/,i), means Di^i < and the term spoils the 
HT scahng equation (10.31). Since the coefficients are rational functions of c, they 
cannot vanish suddenly when c becomes larger than Ccrit.{^,i)- Consequently whenever 
< Ccrit.i^ji) < 2, the coefficient of the corresponding term is expected to vanish 
identically. 

We have checked that this argument is consistent with our previous explicit 
calculations. We have used equations (10.28), (10.8) and the already calculated 
coefficients in equation (10.14) to calculate 

-yi+ Yl - , (10.33) 

i:Xi(i)>i 
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up to order 7. For all the 50 terms with < Ccru. < 2, the corresponding are 
identically zero. 

The existence of an infinite volume limit implies that the small denominator 
problem can be evaded for any c such that < c < 2. We have constructed the 
in terms of the ai. However we could have proceeded directly, writing the a^_^_^^i in 
terms of the a^ f 

at^. , = M\al^^ + ^'i'<k<,, + ■ ■ ■ (10.34) 



k+q>l 



The coefficients ff'^ and the higher order ones can be obtained by using the expansion 
of equation (10.27) in the logarithm of equation (3.12) and expanding order by order 
in a^. The series does not terminate. The linear transformation is the same as before 
because of and ai only differ by nonlinear terms. Using 



< , = n\K , (10.35) 



we obtain 



si, = . (10.38) 



= \i)Ki + E ^'I'KkK, + ■ ■ ■ (10-36) 

k+q>l 

We then introduce the expansion 

i:Xi(i)>i m 

and obtain 

with Nf, given by a formula similar to equation (10.17), except that it does not 
terminate. A detailed analysis shows that the two formulas have in common that the 
numerator depends only on coefficients of strictly lower orders in /3, and equation (10.38) 
can be used order by order in (3 to construct the ; for generic values of c. 

Since TZ~^ is upper triangular, we see from equation (10.35) that hi is equal to ai 
plus terms which go to zero faster. Consequently, for large n, the leading scaling is 

Ki oc A^o . (10.39) 

Following reasonings used before, this implies that terms in the expansion equation 
(10.37) that scale faster than yi for any < c < 2 should have a vanishing coefficient. 
In other words: 

< Cerit. (/, i) <2^s^,i = 0. 

Given the specific form of the s^j given in equation (10.38), the hi have no poles for 
< c < 2. The ai being linear combinations of hi and the a; being linear combinations 
of products of ai, we conclude that the expansion of the a; in terms of the scaling 
variables have also no poles for < c < 2, 

Again we see that there exists a unique continuous definition of the scaling variables 
that can be used at particular values of c where the denominator is exactly zero. From a 
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practical point of view, the calculation at fixed c of the s^; is easier than the calculation 
of the because no limit needs to be taken explicitly. The sf^ being rational function 
of c cannot be zero everywhere except at isolated values. Consequently, we can set to 
zero the s^; having Ccrit.Q'i i) < 2 even at values of c where Di i = 0. 

The initial values yo have a very simple interpretation. We know that yn,i is the 
only leading term of a^ i when n becomes large. If at a given < c < 2, a nonlinear 
terms scales exactly like yn,i, then by increasing c slightly (but keeping c < 2), we can 
make this term dominant in contradiction with the existence of the infinite volume limit. 
Consequently, 



lim Ar"< I = lim A^^^n,/ = yo,i ■ (10.40) 

n— >oo n— >oo 

From equation (10.29), we see that 

yo,i = i-Py^x^"'^ ■ (10.41) 

Furthermore [118], we can consider the x*^^'-* as functions of the initial values oq,;. 
If we now replace these initial values by the m-advanced values am,i, we find that 
X^'^''Ho'm,i) — ^rX^^'H^o,;)- 1^ t^is sense, the infinite volume limit quantities x^^'^ can be 
seen as scaling variables. 

10.6. The nonlinear scaling variables of the nontrivial fixed point 

The construction of the nonhnear scaling variables can be repeated verbatim [119] for 
the non trivial fixed point. However, unlike the HT case, all the calculations have to 
be performed numerically. Starting from the basic quadratic map of equation (10.1), 
introducing new coordinates that are zero at the non-trivial fixed point, reexpressing 
these coordinates as linear combinations of the right eigenvectors and fixing the scale 
(for instance, by requiring that the HT fixed point is located at (1, 1, . . .)), we obtain 
a RG transformation for the linear scaling variables (denoted h) that can be written in 
the form: 

~ 1 + APh +A^^/i h ' ^ ^ ^ 

J. T i\rnji,p ^0 "'n,p"'n,q 

As in the HT case, the first subscript refers to the number of iterations and the second 
is the index of the variable. 

As discussed in section 7, for D = 3 there is one and only one eigenvalue larger 
than 1. We can express the linear scaling variables in terms of the nonlinear scaling 
fields (denoted y): 

Kr -J2*r,illyn:m^ (10-43) 

i m 

and proceed as before. One can also find expansions of the scaling fields in terms of the 
linear variables, by setting 

yn,r^J2''<^Il'^n":m, (10.44) 
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The only potential problem comes from small denominators. In [119, 113], a partial 
numerical survey of possible small denominators was done. The worse case found was 
A2 — A4 with two parts in a thousand. 

A major difference with the HT case, is the existence of a relevant direction. 
Consequently, yi plays a very special role as the coordinate along the unstable direction. 
As yj can be expressed in terms of the linear scaling variables which are themselves linear 
combinations of the original coordinates a, we can have 

|/i(a) = (10.45) 

as the equation defining the stable manifold. Similarly, the unstable manifold in the 
a coordinates corresponds to the 1-dimensional trajectory a(|/i, 0, 0, . . .). For practical 
purpose, one can expand the linear scaling variables to large order in yi and to low order 
in a few irrelevant variables. 



10.7. Convergence issues 

Up to now, it has been shown that it seems possible to construct formal expansion of the 
linear scaling variables in terms of the nonlinear scaling variables or vice-versa, for the 
HT fixed point or the nontrivial fixed point. This does not mean that these expansions 
define analytical functions. On the contrary, the existence of multiple fixed points, 
suggests that these expansions have at best a finite radius of convergence. Numerical 
experiments [119] testing the scaling of the nonlinear scaling variables suggest that 
the two expansions have overlapping region of convergence. This will be illustrated in 
subsection 11.3. 



10.8. The scaling variables of the Gaussian fixed point 

The eigenvalues of the parity preserving linearized RG transformation at the Gaussian 
fixed point [30] are 

Xq. = 2c-^ = e^-jiD-^) (10.46) 

for J = 1, 2, . . .. The interpretation is simply the scaling of the parity invariant couplings 
in g2j(f)^-' interactions. For D > A there is only one relevant direction corresponding to 
the mass term. For 2 < D < 4, there are at least two and at most a finite number of 
relevant directions. 

This problem can be reformulated in term of the evolution operator of the harmonic 
oscillator 

H^^.'^^ (10.47) 

during a finite euclidean time t — —it. The correspondence is c = exp(2ci;r) and 
m = /3(c— l)/(2 — c) 'm.h—1 units. The connection with the £ ~ 1 + 5 limit discussed 
in section 9 is a;T = 5{D — 2)/2. 
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There are many zero denominators in integer dimensions, e.g., Ai = for D = 3. 
If the numerators are not zero, one can modify [31] the situation by considering n- 
dependent coefficients. This generates logarithmic corrections which are necessary. 
These can be observed in the large order of the HT temperature expansion in reference 
[120]. Another possibility is to use the idea of dimensional regularization [121] as already 
explained in the HT case. This might help reinterpreting the connection between the 
1/e poles in the D — A — e regularization and the logarithmic divergences in a cutoff 
regularization [122]. 

Practical constructions of the nonlinear scaling variables remain to be developed. As 

the radius of convergence of perturbative series is zero it is not clear that the procedure 
would work as in the two other cases. Modified perturbative methods where a cutoff in 
field space is introduced [85, 88, 91] might work better. 

We are not aware of any explicit construction of the scaling variables for the low- 
temperature fixed point. 

11. Interpolation between fixed points and critical amplitudes 

11.1. Global RG flows 

The critical amplitudes are in general non-universal and their calculation requires that 
we go beyond the linear approximation and interpolate among fixed points. In the rest 
of this section, we consider mostly the flows from the non-trivial flxed point to the HT 
fixed point. 

Given the success of field theoretical methods based on perturbation theory 
[33, 34, 23], it would certainly be desirable to extend the construction of nonlinear 
scaling variables with initial condition near the Gaussian fixed point. As explained 
in the previous section, this construction of the nonlinear scaling variables remains to 
be done. The basic picture is that the Gaussian fixed point is located on the stable 
manifold of the nontrivial fixed point and that it is possible, for D < 4, to use the 
unstable directions of the Gaussian fixed point to reach the nontrivial fixed point. For 
this reason, the Gaussian fixed point is often called the UV fixed point and the nontrivial 
fixed point the IR fixed point. The situation is depicted in figure 5. 

11.2. Critical amplitudes and RG invariants 

In section 10, we have seen that the construction of the coordinates a in terms of the 
nonlinear HT scaling variables provides a formal solution to the problem of the RG 
fiows. However, this is not the end of the story since it is not clear how accurate finite 
order expansions can be. Also the initial values were identified up to a factor {—f3Y /2r\ 
with the infinite volume susceptibilities x^'^^\ the very quantities that we would like to 
compute! 

Indeed, our goal is to compute x^^^^ or equivalently ^o,r for initial conditions near the 
nontrivial fixed point where it is easy to use the other scaling variables y. For simplicity. 
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Figure 5. A qualitative description of a RG flow starting near the Gaussian fixed 
point, passing by the nontrivial fixed point and ending on the HT fixed point. 



we will assume that the initial conditions are exactly on the unstable direction and close 
to the nontrivial fixed point: 

yo,i = u , 
yo,2 = , 

yo,3 = . . . 

More complicated cases are discussed in reference [119]. Under a RG transformation, 
u ^ Xiu while the HT nonlinear scaling variables transform multiplicatively according 
to equation (10.13). Prom equations (7.23) and (10.6), we can rewrite 

This imphes that 

Cr = ynAVnAV''- , (11-2) 

are n-independent or in other words, RG-invariant. These relations suggest that the 
two fixed points are in some approximate sense dual [123] to each others. We can now 
rewrite the basic quantities that we want to calculate as 

yO,r = yoAyOAV'-iyO,!)-'"' = Criu)--^'^ (11.3) 

The constant Cr are not universal. They depend on the normalization choice for y, 
however equation (11.3) is independent of this choice. It is possible [28] to relate u to 
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Figure 6. Ci{u) defined in equation (11-4) versus u. 



P — Pc for a particular model and consequently, calculating Cr provides an estimate of 
the critical amplitudes. The remaining question is: can we calculate C^? 



11.3. Overlapping regions of convergence 

In principle, the RG invariants can be calculated for any n. In practice, low n 
calculation fail because of the low accuracy of the HT expansion and large n calculations 
fail because of the low accuracy of the expansion in the nontrivial scaling variables. The 
accuracy of the two expansions in intermediate regions can be tested empirically by 
monitoring the stability of the estimates of Cr- In figure 6, we have displayed 

C,{u)^y,{h{hw{u))u\ (11.4) 

with the HT scahng variable yi(h) calculated up to order and h.w{u) calculated 
up to order 80 in u. The RG invariance only implies that Ci{Xu) — Ci{u) (we remind 
that A ~ 1.427). A wide plateau is observed in figure 6 for 1 < u < 3, indicating good 
convergence properties in overlapping regions. RG invariance only forces periodicity on 
a log scale, but apparently the log-periodic oscillations are very small. They can be seen 
better in simplified models [28]. Log-periodic oscillation are discussed into more detail 
below in subsection 11.5. 
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11.4. Approximately universal ratios of amplitudes 

Using the first Imax HT nonfinear scafing variables, it is possible to construct Imax — 1 
constants of motion: 

-P'0!(Z^-j5^ (11.5) 

These quantities are RG invariants. We can evaluate them at n = 0. Using equation 
(10.41), we obtain 

= (-l)'^'(^(2))(.-l)(D/2)+. (11-6) 

If the conjecture [124] that {-ly+^x^'^'^^ > is correct, then Gr > 0. 

We now evaluate the universal ratio on the unstable direction of the unstable fixed 
point. We call Gr{u) the corresponding value, we have 

Gr(Xu) = Gr(u) , (11.7) 

Consequently, we have the Fourier expansion: 

Gr{u) = J2 ^r,,u''^ = E A,,e^^'^"^" , (11.8) 
g g 

with 

The function is clearly periodic in ln(M) and we call the oscillations due to the non-zero 
Fourier modes "log-periodic" . The coefficients can be calculated as 

1 dv 

Ar,, = ^ / —U-^'^^Gr{u) . (11.10) 

The oscillatory terms are very small, as noticed in references [27, 28, 123] and we 
have the approximate universal ratios 

Gr{u) ~ Ar,o . (11.11) 

These constants can be estimated using the methods discussed in subsection 11.3. The 
smallness of the nonzero Fourier modes also applies the non- universal function Gi{u) 
discussed in subsection 11.3. 



11.5. More about log-periodic corrections 

The possibility of log-periodic terms were first discussed in references [32, 45]. They were 
identified in the high-temperature expansion [27, 28] of the HM. The amplitudes Aper. are 
however quite small, typically, they affect the 16-th significant digit of the susceptibility 
and it takes a special effort to resolve them numerically. They are amplified [27, 28] 
by estimators of critical exponents such as the extrapolated slope Sm designed [125] to 
remove subleading corrections in estimation of the critical exponents at successive order 
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Figure 7. The analog of Ci in figure 6, for the quadratic map of equation (11.12) 
with ^ = 0.1, versus the natural logarithm of the scaling variable y. 



in the HT expansion. This prevents an accurate determination of 7 from HT expansion 
[27, 28]. 

Limit cycles often appear in two dimensional ordinary differential equations. 
However, their stability in higher dimensions is an issue debated in the dynamical system 
community. For instance, the Landau scenario for the onset of turbulence, based on the 
appearance of limit cycles, is not considered viable. For the case of interest here, the 
log-periodic oscillations reflect the discrete invariance of the original hamiltonian of 
equation (3.1). This symmetry protects the periodicity even in the continuum limit. 
However, this remnant of the discrete structure is very small numerically. 

It is possible to design toy models [28] where the effect is larger for instance, the 
quadratic map 

hn+l^fhn + (l-Ohl, (11-12) 

for small values of ^. The first nonzero Fourier mode is quite visible in the analog of Ci 
defined in equation (11.4) as shown in figure 7 for ^ = 0.1. 




-10 12 3 

ln(y) 
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12. Nontrivial Continuum limits 

12.1. The infinite cutoff limit 

In this section, we apply the general procedure outlined by Wilson in reference [32] for 
the approximate recursion formula of equation (5.1) to the HM . We consider a sequence 
K = 1,2... of models with (3 = f3c — Xi^u where u is positive but not too large and 
Ai the only relevant eigenvalue as in section 11. (3c depends on the particular choice of 
initial local measure Wq . We introduce the increasing sequence of UV cutoffs 

A;, = 2T>Aj, , (12.1) 
with Ar a scale of reference. We define the renormalized mass 



u 



where x^^^ (/?c — ^u) means the susceptibility a,t P — Pc — ^u. Given that 

A7 = 2*, (12.3) 
the dependence on the UV cutoff disappears at leading order and one obtains 

m% = 7^ (12.4) 

with the log-periodic corrections 

LPC = 41 cos [uj [\nu + ^In (^^^ + 0(2)) . (12.5) 

In the infinite cut-off limit {K oo), the sublcading corrections disappear. On the 
other hand, the LPC do not, and we are in presence of a limit cycle with a cutoff 
dependence quite similar to references [126, 127]. Strictly speaking the infinite cutoff 
limit does not exist, however, for practical purpose, the effects of the oscillations are so 
small that it introduces uncertainties that are smaller than the accuracy with which we 
establish the universality. 

We could now define renormalized coupling constants in cutoff units using the higher 
order %('') just as we have done in equation (12.2) for the renormalized mass. However, 
it is usually more convenient [128] to use the value of these couplings in units of the 
renormalized mass. We thus consider dimensionless couplings of the form 

[/(^) oc x^'^H/Sc - Ar^«)mf . (12.6) 

The fact that the quantity is dimensionless implies that the UV cutoff dependence 
disappears and we are left with a quantity proportional to the ratios of susceptibilities 
defined in equation (11.6). For D = 3, we define 

C/(20 = lim;3_;3e(-l)'+'X^''^ (;^(2))(3-5/)/2 ^3(1-0/2 _ (i2.7) 

In subsection 11.4, it has been argued that along the unstable manifold, this quantity is 
approximately universal. Indeed by taking the limit (5 — > /3c, the flow starts on the stable 
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2/ 


f;(2/)* 


4 


1.505871 


6 


18.10722 


8 


579.970 


10 


35653.8 


12 


3.57769 10^ 


14 


5.3176310^ 


16 


1.0972010^^ 


18 


3.00025 10^^ 


20 


1.0199810^^ 



manifold but then ends up on the unstable manifold and the ^7^^') should be universal 
in the approximation where the very small log-periodic oscillations are neglected. This 
approximate universality means that once we have picked the renormalized mass, all 
the other renormalized couplings are completely fixed. 

A counterpart of this discussion for the LPA of ERGE can be found in sections 
2.10 and 3.4 of reference [8] where a discussion of the various continuum limits that 
can be constructed near the Gaussian fixed point can also be found. Field theoretical 
approaches of IR stable trajectories are also discussed in references [24, 25, 26]. It would 
be interesting to see how the notions developed in these articles (for instance the "large 
river effect") can be used for the HM. 

12.2. Numerical estimates of the universal ratios 

These expectations have been checked numerically in reference [129], by calculating the 
for four different measures. The results were consistent with universality with six 
or seven significant digits. The results are given in table 3 with uncertainties of order 
one in the last printed digit. Various fits of the asymptotic behavior were performed in 
reference [129] and it was concluded that the leading growth is consistent with 

^ q\ . (12.8) 

This is illustrated in figure 8 where ln(t/*^^^"'"^)*/t/*^^^)*) is plotted versus ln(2g). A 
factorial growth as in equation (12.8) would imply a straight line with slope 2, which is 
very close to the slope 2.1 of the hnear fit in figure 8. 

This factorial growth is similar to what is found in reference [130, 131, 132, 133] 
for other models studied in the context of multiparticle production. Note that the 
generating function of the connected 2/-points function has a 1/(2/)! factor at order 21 
(see equation (2.8)) which means that the expansion of the generating function of the 
connected functions in powers of an external field has a finite radius of convergence. 
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1.5 2.0 2.5 3.0 
ln(2q) 

Figure 8. ln([/(29+2)*/{/(29)*) versus ln(2g). The straight line has a slope 2.11. 
12.3. Other universal ratios 

The U"^^* were calculated by picking an arbitrary measure, finding f3c corresponding to 
that measure and then calculating i?„ until the ratios stabilize with sufficient accuracy. 
In this process, it can be noticed that when Rn — R*, typically after 20 or 30 iterations, 
the C/^'" temporarily stabilize at different values than the final ones. We call these 
temporary values U'^'^*. This is illustrated in figure 9 for r = 2. This figure shows the 
importance of going sufficiently far from the nontrivial fixed point to get the correct 
answer. Similar behavior with temporary plateaus can found for the LPA of ERGE in 
reference [26]. U'^'^* can be calculated from the the numerical coefficients of R*. For 
instance, 

U^* ~ 24 X (0.053537 - 0.35871V2) /(2 x 0.35871)^/^ ~ 0.8287 (12.9) 

This number was calculated by expanding the logarithm of equation (7.9) and then 
fixing the normalizations as in equation (12.7). A few other values are given in table 4. 
They are clearly different from the C/^^*. Prom values for r up to 10, their growth seems 
consistent with a simple factorial. 
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Figure 9. U^^^* estimates as the number of iterations n increases. 



Table 4. Universal values [/(^O* _ 



21 




4 


0.828719 


6 


4.17757 


8 


49.3335 


10 


1033.20 



12.4- The critical potential of the symmetric phase 

Following section 2, the effective potential can now be expressed as a function of x^'^\ 
(pc and the C/^''*. To simplify the notations, we define m = (x*-^-')"^''^, the renormalized 
mass in cutoff units. For D = 3, the effective potential reads 

Veff{(j)c) = rrv'F{^,m-"'') . (12.10) 

with 



r=l 



2r 2r 

2r! 



(12.11) 



and /2 = 1, /4 = -U^\ /e = lQ{U^*f - U^'* etc... Equations (7.15) and (12.10) have 
the same form. The correspondence is m ^"""""^ . If we calculate the universal ratios, 
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Table 5. Universal values of r2i calculated numerically and compared to values 
obtained for the nearest neighbor Ising universality class. 



21 


r2i 


r2i for other models 


6 


2.0149752 


2.048(5) [134] 


8 


2.679529 


2.28(8) [134] 


10 


-9.60118 


-13(4) [134] 


12 


10.7681 


20(12) [137] 


14 


763.062 


560(370) [137] 



the parameter m or disappear and we obtain numerical values f/^^')* for equation 

(7.15) and the distinct values for equation (12.10). Consequently, the function F 

is distinct from the function U in equation (7.15). 

In reference [129], rescaling in both coordinates were applied in order to compare 
with a parametrization introduced by Campostrini, Pelissetto, Rossi and Vicari [134] 
where the effective potential is expressed in terms of a universal function 

A{z)^zy2 + z'/24 + Y,^/' , (12.12) 

l>3 ^ 

We have 

T2l = /2.//r' , (12.13) 

and these quantities can be trivially reexpressed in terms of the t/^^'*, for instance, 

77(6)* 

^6 = 10 - pi)^ . (12.14) 

Numerical values in table 5 are not very far from those calculated for nearest neighbor 
models. Approximate relations among the r2i [135, 134] were also checked in reference 
[129]. The first 4 coefficients of the field expansion are positive and there is no convexity 
issue near the origin. The discussion of the convexity at arbitrary field strength, requires 
some understanding of the radius of convergence of the expansion. The convexity of the 
effective potential has been demonstrated for ERGE, using spectral representations for 
the RG flows, in reference [136] where a connection between this issue and the finite 
time singularity discussed above in subsection 9.5 is also drawn. 



13. The large-N limit 

13.1. Calculations at finite N 

If we keep the 0{N) symmetry unbroken, the Fourier transform of the local measure 
depends only on k.k = u. Here A; is a source conjugated to the local field variable 0. 



Nonlinear Aspects of the RG Flows of Dyson's Hierarchical Model 



58 



Replacing khy u and the second derivative by the A^-dimensional Laplacian in equation 
(3.12), we obtain the RG transformation for the Fourier transform of the local measure: 

Rn+iA^) oc cxp ( - l^i^u^ + 2iV^)) (i?„,^(cM/4))' , (13.1) 

The values of the exponents 7, A and the inverse critical temperature for a measure 
generalized Ising measure (or nonlinear sigma model measure) 5(0.0 — 1) calculated in 
reference [18] are given in table 6. To facilitate the comparison, we also display z/ = 7/2 
(since 77 = here) and u; — A/u in table 7. The HM results coincide with the 4 digits 
given in column (2) of table 3 (for u) and 4 (for u) in [101] for Polchinski equation. They 
coincide with the six digits for 1/ given in the line d — 3 of table 8 of [14] for the HM with 
A^= 1, 2, 3, 5 and 10. As in the case = 1 discussed before, we found discrepancies 
of order 1 in the fifth digit of u and slightly larger for uj with the values found in table 
1 of [57]. For = 1, the same discrepancy can be found in references [17, 19]. Our 
estimated errors are of order 1 in the 9-th digit. For A^ = 1, this is confirmed by an 
independent method [56]. For A^ = 2, 3, 5, and 10, this is confirmed up to the sixth digit 
[14]. Consequently, a discrepancy in the 5-th digit cannot be explained by numerical 
errors. It seems clear that the two models are inequivalent. Note also that for A" > 2, a, 
the specific heat exponent shown in table 7, is more negative than for nearest neighbor 
models [23, 34]. 



13.2. Ma's equation 

The basic RG equation in the large A^ hmit was first derived by Ma [138]. It can be used 
for conventional 0{N) sigma models or for the 0{N) version of the HM. We consider 
the partition function 

/■+00 _ _ 

Z{J)=Y[ d^0^e-^+^--^-^% (13.2) 



with 



^ = E ^-^-A + E ^o{<l>l) ■ (13.3) 



xy 



_ — ♦ — ♦ 

We use the notation 0^ = (px-'Px ^-nd A^y is a symmetric matrix with negative 
eigenvalues. We assume that '^^.A^y — 0. This condition is not satisfied by the 
quadratic form appearing in equation (3.1). Using results of subsection 14.1, it is possible 

to show that a term A* 0^ must added to f3H, the non-local part of the action, in 
order to satisfy this condition. In order to keep the original partition function invariant, 
this term must be subtracted from the original potential. Consequently, the potential 
in equation (13.3) is related to the local measure introduced in section 3 by the relation 

H(0') = -lniyo(0) - ■ (13.4) 
One sees that Vq — O for the Gaussian fixed point. Defining the rescaled potential 

Vo{X) = NUoi^) , (13.5) 
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Table 6. 7, A and l3c/N for AT = 1 ... 20. 



N 


7 


A 




1 


1 


29914073 


0.425946859 


1.179030170 


2 


1 


41644996 


0.475380831 


1.236763288 


3 


1 


52227970 


0.532691965 


1.275794011 


4 


1 


60872817 


0.590232008 


1.302790391 


5 


1 


67551051 


0.642369187 


1.322083069 


6 


1 


72617703 


0.686892637 


1.336351901 


7 


1 


76479863 


0.723880426 


1.347244235 


8 


1 


79469274 


0.754352622 


1.355791342 


9 


1 


81827105 


0.779508505 


1.362657559 


10 


1 


83722291 


0.800424484 


1.368284407 


11 


1 


85272636 


0.817977695 


1.372974325 


12 


1 


86561092 


0.832855522 


1.376940318 


13 


1 


87646998 


0.845589221 


1.380336209 


14 


1 


88573562 


0.856588705 


1.383275590 


15 


1 


89372812 


0.866171682 


1.385844022 


16 


1 


90068903 


0.874586271 


1.388107107 


17 


1 


90680338 


0.882027998 


1.390115936 


18 


1 


91221507 


0.888652409 


1.391910870 


19 


1 


91703752 


0.894584429 


1.393524199 


20 


1 


92136121 


0.899925325 


1.394982051 


00 




2 


1 


^ = 1.42366.. 



Table 7. u, oj and a iov N = 1 . . . 20. 



N 


z/ = 7/2 


u = A/u 


a = 2 - 3i/ 


1 


0.649570 


0.655736 


0.051289 


2 


0.708225 


0.671229 


-0.124675 


3 


0.761140 


0.699861 


-0.283420 


4 


0.804364 


0.733787 


-0.413092 


5 


0.837755 


0.766774 


-0.513266 


6 


0.863089 


0.795854 


-0.589266 


7 


0.882399 


0.820355 


-0.647198 


8 


0.897346 


0.840648 


-0.692039 


9 


0.909136 


0.857417 


-0.727407 


10 


0.918611 


0.871342 


-0.755834 
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and using a saddle point calculation of the partition function in the large N limit, it is 
possible to show [138, 139] that = 2dVeff/d(f)l obeys the self-consistent equation 

2C/^(0^ + /a(M2))=M\ (13.6) 

where /a(^^) is the one-loop integral corresponding to the quadratic form A and a 
mass term M^. The prime denotes the derivative with respect to 0^. was denoted 
u in subsection 9.3. For the HM, it is shown in section 8 that the function / 

f"M{z) = Y, , (13.7) 

n=0 \ I ' 

For comparison, for a sharp cutoff model (SCM) in 3 dimensions, we have 

For a generalization of the Ising model (nonlinear sigma models) with a local 
measure 5(0^ — 1) as in reference [18], the saddle point condition simphfies to /(M^) = 
1/N. For the HM, 13c can then be determined from the condition fnui^) = 1/N. 
This imphes Pc — -^(2 — c)/(2(c — 1)). This result was verified in reference [18] and is 
consistent with table 6. 

Let us consider two models, the first one with a rescaled potential Uq, a UV cutoff 
A and a quadratic form A and a second model with a rescaled potential Uq^s, a UV 
cutoff A/S and a quadratic form A5. For D = 3 and in the large- limit, the two 
models have the same dimensionful zero- momentum Green's functions provided that: 

U'oA^') = (13.9) 

^'C/^((0'-/a«(2C/^,5(0^))/5 + /a((2/5^)C/^,5(</''))) 

In the two cases considered above /a = /a<. = / and the fixed point equation becomes 
very simple. In addition, /3 will be set to 1 in the rest of this section. 
Following references [138, 139], we introduce the inverse function: 

F(2f/^((/)2))=02^ (13.10) 

and the function H{z) = F{z) — f{z). With these notations, the fixed point equation 
corresponding to equation (13.9) is simply 

H{z) = SH{z/S^) . (13.11) 

For the SCM, S is allowed to vary continuously in equation (13.11) and the general 
solution is 

F(z)^fscM(z)+Kz^/^ . (13.12) 

For the HM, with D — 3, S can only be an integer power oi £ — 2^/^ and the general 
solution has an infinite number of free parameters: 

F{z) = fHM{z) + J2 K,z'/'+"^^ , (13.13) 
9 
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Figure 10. z versus F^Mi^)- 



with 

a; ^ — ~ 13.6 , (13.14) 

and q runs over positive and negative integers. There exists a unique choice of the Kg 
in equation (13.13) which cancels exactly the singular part of fnM- It was shown in 
reference [140] that the fixed point corresponds to 

1=0 

This expansion has a radius of convergence 2A*c^ ~ 2.7024 for j3 = 1. 
13.3. Singularities of the critical potential 

Following reference [140], we can use equation (13.15) to define F{z) on the negative real 
axis. We remind that F plays the role of cjp'. As wc move toward more negative values 
of F becomes zero within the radius of convergence of the expansion. The situation is 
illustrated in figure 10. Numerically, F*(— 1.5107. . .) = 0. We then reexpand the series 
about that value of z (which corresponds to F = 0^ = 0) and invert it. The resulting 
series is an expansion of 2C/q' in (p. After integration, and up to an arbitrary constant 
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Uq, we obtain a Taylor series for the critical potential U^. We denote the expansion as 



The absolute value of the coefficients appears to grow at an exponential rate. Linear fits 
suggest a radius of convergence of order 2.5. The signs follow the periodic pattern + + 
- -. This suggests singularities along the imaginary axis. This analysis is confirmed by 
an analysis [140] of the poles of Pade approximants. As 0^ exceeds the critical values 
estimated in the previous section, the power series is unable to reproduce the expected 
function Uq. The situation is illustrated in figure 11. For comparison, we have also 
added A*(f)'^ to Uq in order to undo the subtraction of equation (13.4) and obtain a 
function that has a simple relation with the effective potential as in equation (7.15). 
Again, the resulting function appear to be convex. The numerical values of Uq in figure 
11 have been calculated using a parametric representation discussed in reference [140]. 
Finite radius of convergence have been observed in expansions in the power of the fields 
based on ERGE in the LPA [141, 19]. More generally, it would be interesting to compare 
the large- A'" limit for the HM and in the ERGE approach as discussed, for instance, in 
references [142, 143, 144, 101, 145, 146]. 

13.4- Open problems 

The phase diagram of 3D models with a conventional kinetic term shows interesting 
features such as a line of tricritical points ending at the so-called BMB point 
[147, 139, 148]. A similar study should be done for the HM. A general interpretation 
of the complex singularities of the critical potential would be interesting. From the 
numerical values of the exponents, it is possible to estimate the low order coefficients of 
the 1/A^ expansion [149]. The method of calculation of the coefficients using numerical 
values, was developed with the Sterling series for which we were able the first 7 
coefficients accurately. These results are consistent with the hypothesis that the 1/N 
expansions considered are asymptotic but Borel summable (no indication for poles on 
the positive real axis). 

14. The improvement of the hierarchical approximation 

In this section, we describe the HM as a spin model on the 2-adic fine. The main 
motivation for this new way to look at the model is that it suggests a way to modify 
the model in order to approximate nearest neighbor models in D-dimensions. 

14-1. Scalar models on ultrametric spaces 

It is possible to reformulate the HM as a scalar model on the 2-adic line [65]. We 
give here a presentation [66] that does not require a detailed knowledge of the p-adic 
numbers. We will rewrite the hamiltonian of the HM given in equation (3.1) using a 



oo 




(13.16) 



n=0 
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1 2,3 4 5 



Figure 11. UQ{(jP) for the HM with a parametric plot (filled squares), the series 
truncated at order 50 (thick solid line) and Pade approximants [4/1] (thin line slightly 
above the squares) and [5/2] (thin line closer to the squares). The constant has been 
fixed in such a way that the value at the minimum is zero (upper graph). In the lower 
graph, we have added A*^^ to the [5/2] Pade. 



function v{x, y) which specifies the level / at which x and y start to differ. More precisely, 
if X and y are distinct, v{x, y) — I when Xm = ym ior all m such that n > m > I > and 
xi 7^ yi- At coinciding arguments, we define v{x,x) — 0. Referring to figure 1, we can 
see that 2''(^'f) is the size of the smallest block containing both x and y. H can then be 
rewritten as 

H ^ -I{j2^.y<l>.cl>y + Lj2€) , (14.1) 

x,y X 
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where 

K.,= r^4^ ^4^ >^ 4^ (14.2) 

I ii X = y 

and 

L = ((^)-(Jr')(l-^)-^ (14.3) 

As made clear by the above equations, the strength of the interaction between two 
fields (px and (py depends only on the value of v{x,y). Consequently, the invariance of 
v{x, y) under a group of transformation implies the invariance of H under corresponding 
transformations. 

For the reader familiar with the p-adic numbers. The function 2''*^^'^^ is a regularized 
version of the 2-adic distance. Namely, 

I . \ \x — yU ii \x — yU > 1 , , 

2^{^^y) ^ )^ ^'^ ' ^'"^ (14-4) 

1^ 1 if Ix - y|2 < 1 

The 2-adic norm satisfies a relation stronger than the triangle inequality, namely 

k + y|2 < Max(|x|2,|y|2) . (14.5) 

Normed spaces for which this stronger inequality holds are called ultrametric spaces. 
These concepts are explained in more detail for instance in references [63, 61]. 

In order to describe the invariance of v{x, y), we associate with the sequence of O's 
and I's Xn Xi, introduced in section 3 for 2" sites, a rational number of the form 

n 

x = J]a:^2-"^. (14.6) 

m=l 

The reason for this "inversion" is that |2|2 = 2~^ and similarly in momentum space, the 
largest shell will correspond to odd integers, the next shell by the multiple of 2 of these 
odd integers and so on. If two numbers x and y have this form then x + y can also be 
written in this form provide that we drop the integer part of the sum. Equivalently, 
we can write x = g/2" and y = r/2"^ with q and r integers between and 2" — 1 and 
add q and r modulo 2". Since the integers modulo 2" form an additive group, the set 
of fractions associated with the sites form a group for the addition modulo 1. The odd 
integers modulo 2" form a multiplicative group. In the limit n — > oo the group is called 
the 2-adic units. Wc can pick a canonical form for the representatives of such integers 
as 

u = l + 2z (14.7) 

where ^ is a positive integer between and 2"^^ — 1. Obviously, if x has the form of 
equation (14.6) then ux has also this form after discarding its integer part. 

We are now in position to define a group of transformation acting on the fractions 
associated with the sites. If x and a have the form of equation (14.6) and u has the form 
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of equation (14.7), we define a transformation of x depending on a and u and denoted 
a] which reads 

x[u, a] — ux + a , (14.8) 

where the r.h.s is understood modulo 1. It is clear that these transformations form a 
(non-abehan) group that we could call the Poincare group of the HM. We can interpret 
x[0, a] as a translation and x[u, 0] as a rotation like in the usual Poincare groups. In 
that sense, this is a "global" group of transformation. This is in contrast with the 

symmetries noted by Dyson [5] which consists in interchanging Xm xi+ilxi-i xi 

and Xm xi^iOxi^i Xi "locally" which are discussed in section 3. It is possible to 

prove that 

v{x[u,a],y[u,a]) ^v{x,y) . (14.9) 
Since depends only on v{x,y), this implies that 

^xy — '^x[u,a]y[u,a] ■ (14.10) 

that H is invariant under the transformation 

(f>x — 0xKo] (14.11) 
We can then prove that 

< (l>x[u,a](l>y[u,a] > = < 0x02/ > ■ (14.12) 



14-2. Improvement of the hierarchical approximation 

In this subsection, we follow references [20, 21]. We consider a Gaussian model on a finite 
one-dimensional lattice with 2" sites. The Fourier modes of the scalar field are denoted 
where = and k are integers used to express the momenta in |f units. These 
integers are understood modulo 2" in the following (periodicity in momentum space). 
The action reads 

S=^iT.9ik)^'^^-k (14.13) 

k=l 

where g{k) is even, real and positive . 

We proceed in three steps. First, we relabel the momenta in a way which is 
convenient to read their "shell" assignment as in Wilson [3] . This relabeling allows us to 
introduce a group of transformation whose orbits are precisely these shells. The group 
is the multiplicative group of the 2-adic units and its representations are known [150]. 
Completeness can be used to expand the kinetic term , i.e., the function g{k), for each 
of the shells. This solves the bookkeeping problem. Nicely enough, the classification 
of the representations of the group mentioned above comes with an index indicating 
its resolution power (called the degree of ramification). This naturally provides the 
successive orders of our perturbative expansion. We then show that if we only retain 
the trivial representation in the expansion, we obtain the hierarchical approximation. 
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In this limit, the group of transformation is a symmetry of the action which can be 
identified with the symmetry group of the HM mentioned above. 

We can relabel the momenta k. For this purpose, we use a set orthonormal functions 
which is a discrete version of the Walsh system [151]. We first define 

' 1 if A; = -2"-2 + l,....,2"-2-l 

(14.14) 



*o(A;) 



cu^' if A; = -2"-^ 
otherwise 



and 

^i{k) = l-'^o{k) (14.15) 
with the notation cu = 

For a given integer a = Oq + ai2^ + a22^ + + a„_i2"^~^ with ai — or 1 we define 

n-l 

fa{k)^l[^aX2^k) . (14.16) 
It is clear that f*{k) — fa{—k) and we can check that 

k 

A more detailed analysis shows that fa{k) is non-zero only when k — ±k[a\ for a function 
k[a\ which will be specified. More precisely, it is possible to write 

fa{k) = uSk,k[a] + OJ*Sk-k[a] ■ (14.18) 

This relation defines a one-to-one map k[a]. In the following, a will also be treated as 
an integer modulo 2". We can now expand 

2" 

$fe = X!c„/„(A;) . (14.19) 

a=l 

We can then rewrite 

a=l 

where g{a) = g{k[a\). By construction, the Ca are real field variables. For convenience 
we shall also use their complex form 
1 i 

aa=-iCa + C^a) + -iCa-C^a) . (14.21) 

We can now explain the correspondence between this relabeling and Wilson's cell 
decomposition. Clearly, if oq = 1, fa{k) is supported in the high momentum region. 
More precisely, the 0-th shell, i.e, the one integrated first in the RG procedure, consists 
in configurations which can be expanded in terms of the /i+aj2+...(^)- Similarly, the 
modes corresponding tho the l-th shell are made out of the the /2'+ai+i2'+i+....(^)- 
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14-3. The Hierarchical Approximation and its Systematic Improvement 

In the previous section, we have introduced new field variables Ca corresponding to the 
l-th shell when a can be divided by 2' but not by 2'+^. This property is not affected if 
a is multiplied (modulo 2") by any odd number. As explained in subsection 14.1, odd 
numbers form an abelian group with respect to the multiplication modulo 2"^. The orbit 
of this group within the integers modulo 2"' are precisely the sets of numbers that we 
have put in correspondence with the shells. The representations of this group have been 
studied and classified [150]. In order to use this classification, we have to embed the 
labels introduced above and denoted a, in the 2-adic integers. When a can be divided 
by 2' but not by 2'+^, we say that the 2-adic norm, noted |a|2, is 2~^ In the infinite 
volume limit, or in other words when n tends to infinity, the multiplicative group of the 
odd numbers is the 2-adic units. The representations of this group will be denoted Ug. 
This means that if Ui and U2 are 2-adic units, then 115(1^1^2) = Ils{ui)Ils{u2) ■ The label 
s specifies the representation in a way which will be explained below. 

It is easy to construct explicitly the representations Us- A 2-adic unit can be written 
[152] as -u = ±Exp{Az) where z is a (2-adic) integer and Exp the 2-adic exponential. 
Ils{U') is even or odd under multiplication by —1. On the other hand, z is an additive 
parametrization and the z dependence of II^ will be of the form e*^^ where q is an odd 
integer and r a positive integer. Taibleson [150] calls r + 2 the degree of ramification. In 
summary, the label s is a short notation for the parity, r and q an odd integer modulo 



We can now use these representations to expand the the kinetic term function g(a) 
in each shell. For a given shell I, the a have the form 2^ii (so \a\2 — 2~') and we can 
write. 



At finite volume, i.e., at finite n, the units are understood modulo 2"~' and consequently 
the sum over the representations s is restricted to r < n—l—2. The numerical coefficients 
g^^ are easily calculable using the orthogonality relations among the representation. 

The hierarchical approximation is obtained by retaining only the trivial 
representation in the expansion equation (14.22). In this approximation, and using 
the definition introduced in equation (14.21), the action reads 



After a Fourier transform, we obtain a hierarchical model having the general form (with 
arbitrary 6/ couplings at level /) (see section 6). 

The classification of the representations of the 2-adic units suggest that we 
improve the hierarchical approximation by taking into account the additional terms 
in the equation (14.22) order by order in the degree of ramification. Intuitively, this 
corresponds to the fact that the degree of ramification measures the "power of 



2r 




(14.22) 




n-l 



(14.23) 
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resolution" of the representation. Numerically, this works reasonably well: in a simple 
example, the coefficients become smaller as the degree increases [20]. 

The actions written above take a more familiar form after Fourier transformation 
in a. In general, the l-th momentum shell is responsible for interactions among the 
averages of the Fourier transformed fields inside boxes of size 2'. In the hierarchical 
approximation, the interactions depend only on the position of these boxes inside boxes 
of size 2'+^. When the corrections are introduced up to r = Tmax, the interactions 
depend on the position inside boxes of size 2'+^+'"'"'"^ . 

14 -4- The Improvement of the Hierarchical Approximation As a Symmetry Breaking 
Problem 

It is important to realize that in the hierarchical approximation, S is invariant under 
the transformation 

(Ta (Tua (14-24) 

for any odd number u. When the other terms of the expansion are incorporated, this 
symmetry is broken by each term in a definite way. This allows us to use Ward identities 
techniques. We discuss the simplest case below. 

Suppose we want to calculate the 2-point function, using the perturbative expansion 
described in the previous section. First we use the new variables Ca and the inverse of 
the map k[a] defined in section 2 to write 

< <I>fe$„fc >=< c2[,] > . (14.25) 

In the hierarchical approximation, the value of this expression depends only on the 
momentum shell specified by \a\2. In other words, 

< ci >o=< cl >o (14.26) 

where < ... >o means that the quantity is evaluated, at order zero, or in other words, 
in the hierarchical approximation. 

Suppose that we now include a correction 6g{a) to the approximation of g{a). Then 
in first order in this perturbation, we recover the a momentum dependence within the 
shells given by 

< cL cl >, J2(^~g{ub) - 5m) < clcl >o (14.27) 

b 

In the model considered here, the < c^c^ >o contribution can be evaluated 
straightforwardly and we can check that we recover the first term in the expansion 
of {l/g{ua)) — {l/g{a))). The important point is that the corrections are evaluated 
using the unperturbed action. It is clear that similar methods can be used for higher 
point functions and in the interacting case. 

An extension to D-dimension can be constructed easily by noticing the approximate 
correspondence between the integration over successive shells and the block-spin method. 
On a D-dimensional cubic lattice, we can decompose the block-spin procedure into D 
steps (one in each directions). 
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14-5. Other applications 

The 2-adic formulation of the HM can be used for other purposes. For instance, it is 
possible to understand the absence of certain diagrams in the approximate recursion 
formula [4] . The lack of diagrams with three hues having large momenta coming out of 
a vertex can be understood from the fact that the sum of two 2-adic integers with 2-adic 
norm 1 is a 2-adic integer with a strictly lower norm. In simpler words, the sum of two 
odd number is even. 2-adic analysis was also used to study the Symanzik representation 
of the HM [153]. The basic ingredient [64] being that the quadratic form in H defines 
a random walk with a Hausdorff dimension 2/D. 



15. Models with approximate supersymmetry 

In subsections 7.6 and 12.1, we made clear that the continuum limit requires a fine- 
tuning and that this procedure maybe be seen as unnatural. This feature is clearly 
related to the existence of an unstable direction. In four dimensions, there arc several 
ways to get rid of the unstable directions and limit the flow to marginal directions. 
One possibility is to impose a gauge symmetry that forbids a mass term for the gauge 
fields. Another possibility is to introduce new degrees of freedom with opposite statistics 
that partially cancel the quantum fluctuations. This possibility has been exploited in 
the perturbative treatment of supersymmetric models. In the following, we follow this 
second idea and construct models with approximate supersymmetry. Other hierarchical 
models involving fermions can be found in the literature [154, 155]. 

We consider a free action for N massless scalar fields (px^ and fermion fields ipx'^ 
and ipx^: 

Sfree = ^ J] ^^^^^^^^ + E ^^^^^^^^ ' (l^-l) 
x,y,i x,y,i 

where x and y run over the sites and i from 1 to N. The ■^i*^ and ^jJx^ are Grassmann 
numbers integrated with a measure 

n^^x^*^ . (15.2) 

x,i 

We require that D^^ has positive eigenvalues and that we can write 

z 

The free action S^^^'^ -|- Sp*^^ is invariant at first order under the transformation 

bi^=^^Dxy^ (15.4) 



xyH'y 
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The e and e are Grassmann numbers. Integration by part or Leibnitz's rule cannot be 
used for D^y and the order ee variations do not cancel. 



We now give the explicit form for the bosonic part at finite volume. 



sr=-fE(T)" E ( E ^!L„,.....)= (15.5) 

.^BCb_ V- .2 

2(2 — Cb) ^ ^^(.^nmax^-^l)' ' 

with Cb — c = 2^~2/^. The fermionic part reads: 

'^max 

sr=-frE(T)" ( E <„„,.....))( E <L„...-..)) (15.6) 

with Ci? = 2^~^/^. We have the simple relation cb/2 = {cp/^Y. Using the techniques 
explained in [66], one can show that the fermionic operator is the square root of the 
bosonic operator (see equation (15.3)) provided that 

^ = (|^)* (15.8) 

/ — cp / — Cb 

We can introduce local interactions. The Grassmann nature of the fermionic fields 
restricts severely the type of interactions allowed. For instance, for one flavor (N = 1), 
the most general bosonic local measure is 

>V(0, ip, ip) = W{(j)) + ijijA{(j)) (15.9) 

For convenience, we can reabsorb the local quadratic terms, in the local measure. In the 
following, W{4>) will take the Landau- Ginzburg (LG) form: 

W(<l>) oc exp ( - (^1^ + lrnl)<l>' - Xb<P') • (15.10) 

If the two functions W and A arc proportional, the fermionic degrees of freedom 
decouple. The renormalization group transformation takes the form 

W^2Ai.W (15.11) 

A ^2PpAi.W+{—)Wi.W (15.12) 

Cp 

where the * operation is defined as 

{a*b)w . e*(« / rf^M((^?4^)B((*?4^) ■ (15-13) 

The introduction of a Yukawa coupling can be achieved by having A{(f)) linear. Such a 
term breaks explicitly the Z2 symmetry of the LG measure. Models with two fiavors 
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{i — 1, 2) with the type of bihnear couphng appearing in the Wess-Zumino[156] model 
can be written as 

This is not the most general measure, however it closes under the renormalization group 
transformation which takes the form 

W ^ {Wi.T + Ai.A + Bi.B*) = W' 

A ^ (3fW' + —Ai<T 
Cf 

4 

B —B-kT (15.15) 

Cf 

T ^^TicT + (3F—AicT+WF?W'. 
Cf Cf 

In addition, if we impose that the function B has the following form: 

while W, A and T are 0(2)-invariant, the model is then invariant under the R-symmetry 

(0(i) + ,0(2))^e^fl(^{i)+,0{2)) 

iP^^^ e-'^^^^ (15.17) 

We summarize three numerical calculations [157] performed for the second model 
with D — 4. First, the case where the fermions decouple from the bosons was considered. 
W takes the form 

W{<^) oc exp (-((( E('^'^)'+^^(E('^'^)')')) -(IS-IS) 

The value of m|j, defined as the inverse of the zero-momentum two-point function, is 
shown in the top part of figure 12 as a function of m^. These quantities are expressed 
in cut-off units. For reference we have also displayed the one-loop perturbative result 
and the trivial Gaussian result. One sees that the scalar self-interaction moves m|. up 
and m|. ~ 0.2 when goes to zero. The one- loop result is quite good when mj^ is 
large enough but deteriorates when this quantity becomes smaller. 

The second calculation was done in a bosonic model with a bare mass rriB and 
Xb — coupled to a fermion with the following couplings: 

A = {-l-mB)W 

P = QyW (15.19) 

r = ((-l-m^)2 + ^J((0«)2 + (0(2))2))W^. 

The results are shown in the middle part of figure 12 for Qy = VOM ~ 0.28. One sees 
that the Yukawa coupling moves down. For ~ 0.094, mR becomes and for 
smaller of m^, we enter the broken symmetry phase. 
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Figure 12. The renormalizcd mass as a function of the bare mass in a bosonic 0(2) 
model with bare quartic coupling fixed to 0.01 (top), with bare quartic coupling fixed 
to and a Yukawa coupling equal to VO.OS (middle) and with bare quartic coupling 
fixed to 0.01 and a Yukawa coupling equal to \/0.08 (bottom). 
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Finally, the previous calculation was repeated with = 0.01 instead of 0. In 
perturbation theory, the one- loop quadratic divergence cancel when tub — and 



which justifies the choice of coupling constant. The results are shown in figure 12. One 
sees that the Yukawa coupling in part cancels the effects of the scalar self-interaction, 
however, the cancellation is not as good as in the one-loop formula where thr goes to 
zero when mjj goes to zero. Instead, we found numerically that ~ 0.044 when 
goes to zero. It is of course possible to fine-tune Qy in order to get tur — 0. 

16. Conclusions 

In conclusion, we have shown that the calculation of the effective potential at the 
nontrivial fixed point (equation (7.15)) and for a massive theory in the symmetric 
phase (equation (12.10) can be performed numerically with great accuracy. The second 
calculation requires the ability to follow the RG flows all the way from the nontrivial 
point to the HT attractive fixed point. For the HM, this can be accomplished numerically 
or by constructing the nonlinear scaling variables. The situation could be compared to 
a quantum mechanical problem which can be solved accurately and consistently by 
different methods but for which there is no closed form solution. There remain many 
open problems: the construction of the scaling variables near the Gaussian fixed point 
and in the broken symmetry phase, the exploration of the phase diagram for models with 
N components, the explicit calculation of the 1/A^ expansion of the critical exponents 
and most importantly, the improvement of the hierarchical approximation. 

The HM is a good laboratory to explore new ideas. In particular the idea of 
introducing a large field cutoff [85] in order to produce a perturbative series with a 
better large order behavior or to generate a new type of RG ffows when the field cutoff 
is lowered [158] according to the simple scheme 



Our review emphasizes the connection with the ERGE and the possibility to move 
among improvements of the LPA by taking the limit £ — > 1 or discretizing continuous 
RG flow equations. It should be reminded that before approximations are made, ERGE 
can be formulated several equivalent ways. For instance the exact Polchinski equation is 
related to flow equations for the effective action [159, 160, 100] by a Legendre transform. 
However, their derivative expansions are inequivalent and there is room for optimization 
[93]. In contrast, for the HM, despite a program of systematic improvement, exact 
equations for lattice models remain to be found. We hope that some communication 
can be established between these approaches in the future. 

Except for section 15, most of numerical work was done for D = 3. Considering 
models with bosons and fermions producing effects of opposite signs on the effective 



<y 1 



(15.20) 




(16.1) 
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potential may help resolve controversial issues [161] regarding the triviality and the 
stability of the effective potential in D = 4. 

Finally, it would be desirable to apply similar methods for gauge theories. In 
quantum chromodynamics, understanding how asymptotic freedom and confinement can 
be smoothly connected in a single theory amounts to constructing the renormalization 
group (RG) flows of the theory far away from the two flxed points of interest. In this 
example it is expected that nothing dramatic takes place as we interpolate between the 
two regimes, however, understanding confinement in terms of the weak coupling variables 
remains a challenge. Proving the existence of a mass gap in Yang- Mills theories remains 
one of the Clay Millennium Prize Problems. Recent efforts have made to understand 
the basic problems by using decimation procedures on the lattice [38, 40, 41] by using 
ERGE functional methods [37, 39, 10, 41, 42]. As the Large Hadron Collider is almost 
in operation, a more solid understanding of gauge theories should be a priority in the 
theory community. 
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